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Numerical  ocean  modeling  usually  composes  various  initial-  and  boundary-value 
problems.  It  integrates  hydrodynamic  and  thermodynamic  equations  numerically  with 
atmospheric  forcing  and  boundary  conditions  (lateral  and  vertical)  from  initial  states  of 
temperature,  salinity  and  velocity.  Past  observations,  historical  datasets  and 
climatological  datasets  of  the  ocean  have  contributed  greatly  to  the  knowledge  of  the  data 
fields  of  initial  condition,  atmospheric  forcing  and  boundary  conditions.  Change  in  either 
initial  or  boundary  condition  leads  to  a  variety  of  model  solutions.  It  is  necessary  to 
specify  realistic  initial  and  boundary  conditions  to  achieve  better  understanding  and 
prediction  of  the  ocean  behavior.  However,  uncertainty  often  exists  in  both  initial  and 
boundary  conditions.  Up  to  now,  most  studies  on  ocean  predictability  have  usually  been 
for  one  particular  type  of  model  input  uncertainty  within  the  three  types  of  uncertainty 
(initial  conditions,  open  boundary  conditions,  atmospheric  forcing  function).  This  thesis 
investigates  the  response  of  ocean  model  to  the  three  types  of  model  input  uncertainty 
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I.  INTRODUCTION 


Numerical  ocean  modeling  usually  composes  various  initial-  and  boundary-value 
problems.  It  integrates  hydrodynamic  and  thermodynamic  equations  numerically  with 
atmospheric  forcing  and  boundary  conditions  (lateral  and  vertical)  from  initial  states  of 
temperature  (T),  salinity  (S)  and  velocity.  Past  observations,  historical  datasets  and 
climatological  datasets  of  the  ocean  have  contributed  greatly  to  the  knowledge  of  the  data 
fields  of  initial  condition,  atmospheric  forcing  and  boundary  conditions.  Change  in  either 
initial  or  boundary  condition  leads  to  a  variation  of  model  solutions  (Lorenz,  1963;  Chu, 
1999b).  It  is  necessary  to  specify  realistic  initial  and  boundary  conditions  to  achieve 
better  understanding  and  prediction  of  the  ocean  behavior.  However,  uncertainty  often 
exists  in  both  initial  and  boundary  conditions. 

The  first  difficult  problem  in  ocean  modeling  (regional  and  basin-scale)  is  that 
initial  velocity  field  is  usually  not  available  due  to  an  insufficient  number  of  velocity 
observations.  A  widely  used  model  initialization  is  the  diagnostic  mode,  which  integrates 
the  model  from  known  T,  S,  such  as  climatological  data  (Tc,  Sc)  and  zero  velocity  fields, 
while  holding  Tc  and  Sc  unchanged.  After  a  period  (about  30  days)  of  the  diagnostic  run, 
the  velocity  field  (Vc)  is  established,  and  Tc,  Sc  and  Vc  fields  are  treated  as  the  initial 
conditions  for  numerical  prognostic  modeling.  Recently,  Chu  and  Lan  (2003)  found  that 
during  the  diagnostic  initialization  period,  unrealistic  heat  and  salt  ‘source/sink’  tenns  are 
generated  at  each  time  step. 

The  second  difficult  problem  in  regional  ocean  modeling  is  the  uncertainty  of  the 
open  boundary  condition  (Chu  et.  ah,  1997).  At  open  boundaries  where  the  numerical 
grid  ends,  the  fluid  motion  should  be  unrestricted  since  ideal  open  boundaries  are 
transparent  to  motions.  Two  approaches,  local-type  and  inverse-type,  are  available  for 
detennining  open  boundary  condition  (OBC).  The  local-type  approach  determines  the 
OBC  from  the  solution  of  the  governing  equations  near  the  boundary.  The  problem  now 
becomes  selecting  from  a  set  of  ad  hoc  OBCs.  Since  any  ad  hoc  OBC  will  introduce 
inaccuracies  into  a  numerical  solution  (Chapman,  1985),  it  is  important  to  choose  the  best 
one  from  ad  hoc  OBCs  for  a  particular  ocean  model.  Using  a  barotropic  coastal  ocean 
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model  Chapman  (1985)  evaluated  several  of  the  most  used  ad  hoc  OBCs  (clamped, 
sponge,  radiation)  and  found  that  the  best  OBC  consists  of  a  sponge  at  the  outer  edge  of 
the  model  domain  with  an  Orlanski  radiation  condition  (Orlanski,  1976)  while 
determining  that  the  clamped  OBC  is  probably  the  worst  choice.  Applying  these  results  to 
other  ocean  models  needs  further  investigation.  The  local  approach  suffers  drawbacks 
that  may  restrict  its  use:  no  observational  data  considered  and  the  ill-posedness  of  the 
primitive  equations  model  with  ad  hoc  OBC;  in  other  words,  it  is  hard  to  prove  the 
existence  of  a  unique  solution  (Bennett,  1992;  Oliger  and  Sundstrom,  1978).  To  improve 
the  local  approach  by  using  observations  at  open  boundaries,  Shulman  and  Lewis  (1995) 
proposed  a  method  for  determining  OBCs  of  the  shallow  water  model.  Their  method  is 
based  on  the  integration  of  governing  equations  forward  in  time  and  the  selection  of 
OBCs  via  a  specific  inverse  problem  that  minimizes  a  measure  of  difference  (energy 
flux)  between  the  values  of  observed  and  predicted  variables  at  open  boundaries.  Thus, 
their  method  helps  in  selecting  proper  ad  hoc  OBC  by  using  observations  at  the  open 
boundaries. 

Without  any  ad  hoc  OBC,  the  inverse-type  approach  can  determine  the  OBC  from 
the  "best"  fit  between  model  solutions  and  interior  observations.  The  most  popular 
scheme  for  this  approach  is  an  adjoint  method,  which  consists  of  four  elements:  a  set  of 
control  parameters  or  a  control  vector  (e.g.,  the  unknown  OBC),  a  numerical  ocean 
model,  a  cost  function,  and  an  adjoint  equation.  The  cost  function  is  usually  defined  by 
the  difference  between  observations  and  their  model  counterparts.  The  adjoint  equation  is 
derived  from  minimizing  the  cost  function  with  respect  to  the  control  parameters.  The 
advantage  of  using  the  adjoint  method  is  the  well  posedness  and  the  use  of  observational 
data.  Seiler  (1993)  successfully  detennines  the  unknown  OBCs  for  a  quasi-geostrophic 
ocean  by  using  the  adjoint  method.  The  disadvantages  that  may  restrict  its  use  are  the 
requirement  of  large  amounts  of  computer  time  and  memory,  problems  of  stable 
integration  of  the  adjoint  equation,  the  ocean-model  dependency  of  the  adjoint  equation, 
and  difficulty  in  deriving  the  adjoint  equation  when  the  model  contains  rapidly  changing 
processes,  such  as  ocean  mixed  layer  dynamics.  Chu  et  al.  (1997)  propose  a  simplified 
method  that  overcomes  the  disadvantage  of  the  current  inverse-type  approach.  This 
method  can  determine  OBCs  of  any  ocean  model  (i.e.,  a  universal  method)  from  interior 
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observations.  The  essence  of  the  method  is  to  seek  the  relationship  among  three  vectors: 
open  boundary  parameter  vector  (B),  observation  vector  (O),  and  solution  vector  (S).  If  B 
is  given,  we  can  integrate  the  numerical  ocean  model  and  get  the  solution  vector  (S).  If  B 
is  unknown,  the  optimization  method  is  used  to  detennine  B  by  minimizing  the  root- 
mean  square  difference  between  O  and  S. 

The  third  difficult  problem  is  the  uncertainty  in  the  atmospheric  forcing  function. 
This  is  largely  due  to  the  lack  of  meteorological  observations  over  the  ocean  surface.  For 
example,  Chu  et  al.  (1999c)  found  significant  difference  in  wind  forcing  over  the  South 
China  Sea  during  the  lifetime  of  tropical  cyclone  Ernie  (November  4-18,  1996)  between 
NASA’s  Scatterometer  (NSCAT)  and  NOAA’s  National  Centers  for  Environmental 
Prediction  (NCEP)  winds.  The  root-mean-square  difference  increased  from  3.6  m  s'1  on 
November  1  to  a  maximum  value  of  6.7  m  s'1  on  November  4,  1996,  which  was  the  day 
the  boundary  current  was  strongest,  fluctuating  afterward  between  6.7  and  2.7  m  s’1.  It 
varied  from  50  to  100%  of  the  internal  wind  variability  and  equaled  the  internal  wind 
variability  on  November  12.  Chu  et  al.  (1999c)  further  investigated  the  model  uncertainty 
due  to  the  uncertainty  of  the  surface  boundary  conditions  using  the  Princeton  Ocean 
Model  (POM)  with  20  km  horizontal  resolution  and  23  sigma  levels  confonning  to  a 
realistic  bottom  topography  during  the  lifetime  of  tropical  cyclone  Ernie  (November  4- 
18,  1996).  The  uncertainty  in  surface  winds  generated  uncertainty  in  currents  and 
thennohaline  structure.  The  root-mean-square  difference  of  currents  between  NSCAT 
and  NCEP  runs  decreased  with  depth,  increased  with  time,  and  rose  to  a  maximum  value 
of  18  cm  s’1,  occurring  at  the  surface  on  November  30.  This  varied  from  47  to  102%  of 
the  internal  variability  of  the  surface  currents.  The  uncertainty  in  surface  winds  generated 
uncertainty  in  temperature.  The  root-mean-square  difference  of  temperature  between 
NSCAT  and  NCEP  runs  increased  in  depth  from  the  surface  to  a  subsurface  level  (sigma 
=-0.025)  by  a  maximum  value  of  0.52°C. 

Within  the  nonlinear  dynamics  community,  three  types  of  sensitivity  are 

commonly  discussed:  sensitivity  to  initial  conditions,  sensitivity  to  parameters,  and 

sensitivity  to  boundary  conditions.  The  first  type  of  sensitivity  has  received  a  great  deal 

of  attention,  the  second  has  received  a  good  deal  of  attention,  but  the  third  has  received 

relatively  little  attention.  Recently,  Chu  (1999b)  found  that  two  kinds  of  predictability 
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problems  exist  in  the  Lorenz  system,  namely,  the  model  sensitivity  to  initial  (first  kind) 
and  boundary  (second  kind)  perturbations.  The  effect  of  the  boundary  error  on  the  model 
can  be  represented  as  a  forcing  term. 

Introducing  the  same  small  relative  error  (10‘4)  to  either  the  initial  or  boundary 
condition,  the  Lorenz  system  has  a  growing  period  and  an  oscillation  period.  During  the 
growing  period,  the  model  error  increases  from  0  to  an  evident  value  larger  than  1. 
During  the  oscillation  period,  the  model  error  oscillates  between  two  evident  values.  For 
the  wide  range  of  the  parameter  space,  both  the  error  growing  period  and  the  relative 
error  are  comparable  between  the  first  kind  of  problem  (initial  inaccuracy)  and  the  second 
kind  of  problem  (boundary  inaccuracy).  This  suggests  the  importance  of  preparing 
accurate  boundary  conditions  for  numerical  prediction.  The  Lorenz  system  is  a  low  order 
(three  components)  convective  model.  As  a  result  the  following  question  arises:  do  the 
two  types  of  predictability  in  Lorenz  system  also  exit  in  realistic  regional  ocean  models? 
The  boundaries  in  ocean  models  include  upper  boundary  (atmospheric  forcing)  and 
lateral  boundary.  Ocean  model  uncertainty  caused  by  the  errors  in  the  atmospheric 
forcing  functions,  such  as  surface  winds  and  thermohaline  fluxes  and/or  in  the  lateral 
open  boundary  conditions  is  the  second  kind  of  predictability.  In  ocean  models,  the  two 
kinds  of  uncertainty  may  occur  together.  The  effect  of  joint  initial/boundary  condition 
errors  on  the  model  perfonnance  is  defined  as  the  third  kind  of  predictability. 

Up  to  now,  most  studies  on  ocean  predictability  have  usually  been  for  one 
particular  type  of  model  input  uncertainty  within  the  three  types  of  uncertainty  (initial 
conditions,  open  boundary  conditions,  atmospheric  forcing  function).  This  thesis 
investigates  the  three  types  of  model  input  uncertainty  simultaneously  using  Princeton 
Ocean  Model  (POM)  for  the  Japan/East  Sea  (JES). 

It  is  organized  into  the  following  chapters  in  this  thesis.  In  chapter  II,  the 
oceanography  of  the  Japan/East  Sea  (JES)  is  described.  In  Chapter  III,  the  Princeton 
Ocean  Model  (POM)  is  introduced.  In  Chapter  IV,  the  experimental  design  is  described. 
In  Chapters  V,  the  analysis  methods  are  introduced.  In  Chapter  VI,  model  errors  due  to 
input  uncertainty  are  investigated.  Finally,  in  Chapter  VII,  the  conclusions  are  presented. 
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II.  JAPAN/EAST  SEA  (JES)  OCEANOGRAPHY 


The  Japan  Sea,  known  as  the  East  Sea  in  Korea,  has  steep  bottom  topography 
(Fig.  1)  that  makes  it  a  unique  semi-enclosed  ocean  basin  overlaid  by  a  pronounced 
monsoon  surface  wind.  The  Japan/East  Sea,  hereafter  referred  to  as  JES,  covers  an  area 
of  10  km".  It  has  a  maximum  depth  in  excess  of  3,700  m,  and  is  isolated  from  open 
oceans  except  for  small  (narrow  and  shallow)  straits.  The  JES  connects  with  the  North 
Pacific  through  the  Korea/Tsushima  and  Tsugaru  Straits  and  with  the  Okhotsk  Sea 
through  the  Soya  and  Tatar  Straits.  In  addition,  the  JES  contains  three  major  basins  called 
the  Japan  Basin  (JB),  Ulleung/Tsushima  Basin  (UTB),  and  Yamato  Basin  (YB);  it  also 
has  a  high  central  plateau  called  the  Yamato  Rise  (YR).  The  JES  is  of  great  scientific 
interest  as  a  miniature  prototype  ocean.  Its  basin-wide  circulation  pattern,  boundary 
currents,  Subpolar  Front  (SPF),  mesoscale  eddy  activities  and  deepwater  formation  are 
similar  to  those  in  a  large  ocean. 


A.  THE  THERMOHALINE  STRUCTURE 

The  thermohaline  structure  of  JES  has  been  studied  by  many  investigators  (Gong 
and  Park,  1969;  Isoda  and  Saitoh,  1993;  Isoda  et  ah,  1991;  Maizuru  Mar.  Observ.,  1997) 
using  limited  data  sets.  For  example,  after  analyzing  satellite  infrared  (IR)  images  and 
routine  hydrographic  survey  data  (by  the  Korea  Fisheries  Research  and  Development 
Agency)  for  the  western  part  of  the  JES  in  the  winter  and  the  spring  1987,  Isoda  and 
Saitoh  (1993)  found  that  a  small  meander  of  a  thermal  front  originates  from  the 
Korean/Tsushima  Strait  near  the  Korean  coast  gradually  growing  into  an  isolated  warm 
eddy  with  a  horizontal  scale  of  100  km.  The  warm  eddy  moves  slowly  northward  from 
spring  to  summer. 

Although  the  seasonal  thennal  variability  on  150  m  depth  is  weaker  than  on  the 
surface,  SPF  still  occurs  at  around  40°N  consistently  throughout  the  year,  and  it  is  located 
at  almost  the  same  location  as  at  the  surface.  It  divides  the  water  masses  with  different 
characteristics.  North  of  the  SPF,  the  temperature  is  uniformly  cold  (1°-3°C)  throughout 
the  year.  South  of  the  SPF,  the  temperature  changes  from  5°C  to  9°C.  The  SPF 
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meandering  at  131°E,  134°E,  and  138°E  forms  several  mesoscale  eddies  (Chu  et  al., 
200 1  a, b).  The  SPF  meandering  near  Okin  Gunto  (134°  E)  in  spring  was  previously 
reported  by  Isoda  and  Saitoh  (1993). 

With  limited  data,  Miyazaki  (1953)  found  a  low  salinity  layer  in  the  SPF  region. 
Later  on  Kim  and  Chung  (1984)  found  very  similar  property  in  UTB  which  they  called 
the  JES  Intermediate  Water  (JIW).  After  analyzing  the  comprehensive  hydrographic  data 
for  the  whole  JES  collected  by  the  Japan  Meteorological  Agency,  the  Maizuru  Marine 
Observatory,  and  the  Hydrographic  Department  of  the  Japan  Maritime  Safety  Agency, 
Senjyu  (1999)  demonstrates  the  existence  of  a  salinity  minimum  (SMIN)  layer  (i.e.,  JIW) 
between  the  TWC  Water  and  the  JES  Proper  Water.  The  southwestern  JES  west  of  132°E 
is  the  upstream  region  of  JIW.  The  lowest  salinity  and  the  highest  oxygen  concentration 
are  found  in  the  38°-40°N  areas  west  of  132°E.  The  JIW  takes  two  flow  paths:  an 
eastward  flow  along  the  SPF  and  a  southward  flow  parallel  with  the  Korean  coast  in  the 
region  west  of  132°E.  Analyzing  the  hydrographic  collected  from  an  international 
program,  Circulation  Research  of  the  East  Asian  Marginal  Seas  (CREAMS),  Kim  and 
Kim  (1999)  discovered  the  high  salinity  water  with  high  oxygen  in  the  eastern  JB  (i.e., 
north  of  SPF)  naming  it  the  High  Salinity  Intermediate  Water  (HSIW). 

Recently,  Chu  et  al.  (1998;  1999a)  reported  the  seasonal  occurrence  of  JES  eddies 

from  the  composite  analysis  of  the  U.S.  National  Center  for  Environmental  Prediction’s 

(NCEP)  monthly  SST  fields  (1981-1994).  For  example,  they  identified  a  wann  center 

appearing  in  late  spring  in  the  East  Korean  Bay.  Chu  et  al.  (200 1  a, b)  further  reported  the 

seasonal  variation  of  the  thermohaline  structure  and  inverted  circulation  from  the  Navy's 

unclassified  Generalized  Digital  Environmental  Model  (GDEM)  temperature  and  salinity 

data  on  a  0.5°  x  0.5°  grid  using  the  P-vector  method  (Chu,  1995).  The  GDEM  for  the  JES 

was  built  on  136,509  temperature  and  52,572  salinity  (1930-1997)  historical  profiles.  A 

three-dimensional  estimate  of  the  absolute  geostrophic  velocity  field  was  obtained  from 

the  GDEM  temperature  and  salinity  fields  using  the  P-vector  method.  The  climatological 

mean  and  seasonal  variabilities  of  the  thermohaline  structure  and  the  inverted  currents, 

such  as  the  SPF,  the  mid-level  (50  to  200  m)  salty  tongue,  the  Tsushima  Warm  Current 

(TWC)  and  its  bifurcation,  were  identified.  Using  the  data  collected  from  Conductivity- 

Temperature-Depth  (CTD)  and  Acoustic  Doppler  Current  Profilers  (ADCP) 
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measurements  in  the  southwestern  JES  from  March  to  June  1992,  Shin  et  al.  (1995,  1996) 
found  a  dipole  structure  of  gyres  with  an  anticyclonic  eddy  near  the  Korean  coast  and  a 
cyclonic  eddy  in  the  UTB. 


B.  CURRENT  SYSTEMS 

Most  of  the  nearly  homogeneous  water  in  the  deep  part  of  the  basin  is  called  the 
Japan  Sea  Proper  Water  (Moriyasu,  1972)  and  is  of  low  temperature  and  low  salinity. 
Above  the  Proper  Water,  the  Tsushima  Warm  Current  (TWC),  dominating  the  surface 
layer,  flows  in  from  the  East  China  Sea  through  the  Korea/Tsushima  Strait  carrying  warm 
water  from  the  south.  The  Liman  Cold  Current  (LCC)  carries  cool  fresh  surface  water 
from  the  north  and  northeast  (Seung  and  Kim,  1989;  Holloway  et  al.,  1995).  The 
properties  of  this  surface  water  are  generally  believed  to  be  detennined  by  the  strong 
wintertime  cooling  coupled  with  fresh  water  input  from  the  Amur  River  and  the  melting 
sea  ice  in  Tatar  Strait  (Martin  and  Kawase,  1998).  The  LCC  flows  southward  along  the 
Russian  coast,  beginning  at  a  latitude  slightly  north  of  Soya  Strait,  terminating  off 
Vladivostok  (Fig. 2),  and  becoming  the  North  Korean  Cold  Current  (NKCC)  after 
reaching  the  North  Korean  coast  (Yoon,  1982a). 

The  TWC  separates  into  two  branches,  which  flow  through  the  western  and 
eastern  channels  of  the  Korea/Tsushima  Strait  (Kawabe,  1982a, b;  Hase  et  al.,  1999).  The 
flow  through  the  eastern  channel  closely  follows  the  Japanese  Coast;  it  is  called  the 
Nearshore  Branch  (Yoon,  1982b)  or  the  first  branch  of  TWC  (FBTWC)  (Hase  et  al., 
1999).  The  flow  through  the  western  channel  is  called  the  East  Korean  Warm  Current 
(EKWC),  which  closely  follows  the  Korean  coast  until  it  separates  near  37°N  into  two 
sub-branches.  The  western  sub-branch  moves  northward  and  fonns  a  cyclonic  eddy  over 
UTB  off  the  eastern  Korean  coast.  The  eastern  sub-branch  flows  eastward  to  the  western 
coast  of  Hokkaido  Island,  and  becomes  the  second  branch  of  the  TWC  (SBTWC). 

The  NKCC  meets  the  EKWC  at  about  38°N  with  some  seasonal  meridional 
migration.  After  separation  from  the  coast,  the  NKCC  and  the  EKWC  converge  fonning 
a  strong  front  that  stretches  to  the  zonal  direction  across  the  basin.  The  NKCC  makes  a 
cyclonic  recirculation  gyre  in  the  north,  while  most  of  the  EKWC  flows  out  through  the 
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Tsugaru  and  Soya  Straits  (Uda,  1934).  The  fonnation  of  NKCC  and  the  separation  of 
EKWC  are  due  to  local  forcing  by  wind  and  buoyancy  flux  (Seung  and  Nam,  1992). 
Large  meanders  associated  with  warm  and  cool  eddies  develop  along  the  front. 


C.  ATMOSPHERIC  FORCING 

The  Asian  monsoon  strongly  affects  the  thermal  structure  of  the  JES.  During  the 
winter  monsoon  season,  a  very  cold  northwest  wind  blows  over  the  JES  (Fig.  3a)  as  a 
result  of  the  Siberian  High  Pressure  System  with  a  mean  surface  wind  speed  between  10 
and  15  ms'1.  By  late  April,  numerous  frontally  generated  events  occur  making  late  April 
and  May  highly  variable  in  tenns  of  wind  speeds  and  the  amount  of  clouds.  During  this 
period,  storms  originating  in  Mongolia  may  cause  strong,  warm  westerlies  (Fig.  3b).  By 
late  May  and  early  June,  a  summer  surface  atmospheric  low-pressure  system  begins  to 
form  over  Asia.  Initially  this  low-pressure  system  is  centered  north  of  the  Yellow  Sea 
producing  westerly  winds.  In  late  June,  this  low  begins  to  migrate  to  the  west  setting  up 
the  southwest  monsoon  that  dominates  the  summer  months.  The  winds  remain  variable 
through  June  until  the  Manchurian  Low  Pressure  System  strengthens.  Despite  the  very 
active  weather  systems,  the  mean  surface  wind  speed  over  the  JES  in  summer  (Fig.  3c)  is 
between  3  and  4  m/s,  much  weaker  than  in  winter  (Fig.  3a).  By  July,  however,  high 
pressure  (the  Bonin  High)  to  the  south  and  the  low  pressure  over  Manchuria  produce 
southerly  winds  carrying  wann,  moist  air  over  the  East  China  Sea/Y ellow  Sea.  In 
summer,  warm  air  and  strong  downward  net  radiation  stabilize  the  upper  layer  of  the  JES 
causing  the  surface  mixed  layer  to  shoal.  October  (Fig.  3d)  is  the  beginning  of  the 
transition  to  winter  conditions.  The  southerly  winds  weaken  and  the  sea  surface  slope 
establishes  its  winter  pattern. 
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Figure  1 .  The  Japan  /  East  Sea  geography  and  bottom  topography. 
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Figure  2.  Current  system  in  the  Japan/East  Sea  (from  Tomczak  and  Godfrey,  1994). 
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Figure  3.  Climatological  wind  stress  from  the  COADS  data. 
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III.  THE  PRINCETON  OCEAN  MODEL  (POM) 


Coastal  oceans  and  semi-enclosed  seas  are  marked  by  extremely  high  spatial  and 
temporal  variability  that  challenges  the  existing  predictive  capabilities  of  numerical 
simulations.  The  Princeton  Ocean  Model  (POM)  is  a  time-dependent,  primitive  equation 
circulation  model  rendered  on  a  three-dimensional  grid  that  includes  realistic  topography 
and  a  free  surface  (Blumberg  and  Mellor,  1987). 

In  this  study,  the  model  contains  91  x  100  x  23  fixed  grid  points.  The  horizontal 
spacing  is  10'  latitude  and  longitude  (approximately  11.54  to  15.18  km  in  the  zonal 
direction  and  18.53  km  in  the  meridional  direction)  with  23  sigma  levels  in  vertical 
coordinate.  The  model  domain  extends  from  35.0°N  to  51.5°N  and  from  127. 0°E  to 
142. 0°E.  The  bottom  topography  (Fig.  1)  is  obtained  from  the  Naval  Oceanographic 
Office's  Digital  Bathymetry  Data  Base  5'  x  5'  resolution  (DBDB5).  The  horizontal 
friction  and  mixing  are  modeled  using  the  Smagorinsky  (1963)  form  with  the  coefficient 
chosen  to  be  0.2  for  this  application.  The  bottom  stress  Tb  is  assumed  to  follow  a 
quadratic  law: 

n  =  p,CD  \K\Vb  (1) 

where  p o(=  1025  kg/m  )  is  the  characteristic  density  of  the  sea  water,  V/,  is  the  horizontal 
component  of  the  bottom  velocity,  and  Cd  is  the  drag  coefficient  specified  as  0.0025 
(Blumberg  and  Mellor,  1987).  This  coefficient  is  similar  to  the  value  (0.002)  as  used  by 
Hogan  and  Hurlburt  (2000a, b).  The  significant  attributes  of  the  model  are  as  follows: 


A.  SIGMA  COORDINATE 

It  is  a  sigma  coordinate  model  in  which  the  vertical  coordinate  is  scaled  on  the 
water  column  depth.  The  sigma  coordinate  equations  are  based  on  the  following 
transformation: 


x*  =  x,  y*  =  y,  <j 


Z  -  77  * 

- — ,  t*  =  t 

H  +  fj 


(2) 
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where  x,y,z  are  the  conventional  cartesian  coordinates;  D  =  H  +  77  where  H  (x,  y  )  is  the 
bottom  topography  and  rj(x,  y,  t )  is  the  surface  elevation.  Thus,  cr  ranges  from  a  =  0  at 
z  =  r]  to  <j  =  -1  at  z  =  -H.  The  sigma  coordinate  system  is  necessary  in  dealing  with 
significant  topographical  variability,  such  as  that  encountered  over  continental  shelf 
breaks  and  slopes.  Together  with  the  turbulence  sub-model,  the  model  produces  realistic 
bottom  boundary  layers,  which  are  important  in  coastal  waters  (Blumberg  and  Mellor, 
1987). 


B.  SURFACE  FORCING  FUNCTIONS 


The  atmospheric  forcing  for  the  JES  application  of  POM  includes  mechanical  and 
thennohaline  forcing.  The  wind  forcing  is  depicted  by 


P0K 


M 


du  dv 
dz  ’  dz 


:=0  ~  (T0x,T0y) 


(3) 


where  KM  is  the  vertical  mixing  coefficient  for  momentum,  (u,  v)  and  (rox,  %)  are  the  two 
components  of  the  water  velocity  and  wind  stress  vectors,  respectively.  The  wind  stress  at 
each  time  step  is  interpolated  from  monthly  mean  climatological  wind  stress  from 
Comprehensive  Ocean- Atmosphere  Date  Set  (COADS  1945-1989),  with  a  resolution  of 
1°  x  1°.  The  COADS  wind  stress  was  interpolated  into  the  model  grid  with  a  resolution 
of  10'. 


Surface  thermal  forcing  is  depicted  by 


kJ-1 

H  dz 


=  oc, 


'  Qh^ 


pQ 


+  cc  2  C  (0 0BS  0) 


pj 


(4) 


Ks  —  =  —ax  FS  +  a2  C  (S0BS  -  S)  (5) 

dz 

where  KH  and  Ks  are  the  vertical  mixing  coefficients  for  heat  and  salt,  ( 6 ,  S)  and  (6obs, 
Sobs )  are  the  modeled  and  observed  potential  temperature  and  salinity,  and  cp  is  the 
specific  heat.  The  relaxation  coefficient  C  is  the  reciprocal  of  the  restoring  time  period 
for  a  unit  volume  of  water.  The  parameters  (ai,  0.2)  are  (0,  1)  switches:  ai  =  1,  012  =  0, 
would  specify  only  flux  forcing  is  applied.  Flux  forcing  is  used  in  this  study. 
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C.  LATERAL  BOUNDARY  FORCING 

Boundary  conditions  for  closed  lateral  boundaries,  i.e.,  the  modeled  ocean 
bordered  by  land,  were  defined  using  a  free-slip  condition  for  velocity  and  a  zero  gradient 
condition  for  temperature  and  salinity.  Thus,  no  advective  or  diffusive  heat,  salt  or 
velocity  fluxes  occur  through  these  boundaries.  At  open  boundaries,  the  numerical  grid 
tenninates,  but  the  fluid  motion  is  unrestricted.  Uncertainty  at  open  boundaries  makes 
marginal  sea  modeling  difficult. 


Month 

Feb. 

Apr. 

Jun. 

Aug. 

Oct. 

Dec. 

Tatar  strait  (inflow) 

0.05 

0.05 

0.05 

0.05 

0.05 

0.05 

Soya  strait  (outflow) 

-0.1 

-0.1 

-0.4 

-0.6 

-0.7 

-0.4 

Tsugaru  strait  (outflow) 

-0.25 

-0.35 

-0.85 

-1.45 

-1.55 

-1.05 

Tsushima  strait  (inflow) 

0.3 

0.4 

1.2 

2.0 

2.2 

1.4 

Table  1  The  bi-monthly  variation  of  volume  transport  (unit:  Sv,  1  Sv  =  106  m3  s'1). 

The  bi-monthly  volume  transports  through  open  boundaries  specified  from 
historical  data  are  listed  in  Table  1.  Positive  (negative)  values  are  referred  to  inflow 
(outflow).  Warm  water  enters  the  JES  through  the  Korea/Tsushima  Strait  with  the  TWC 
from  the  East  China  Sea  and  exits  the  JES  through  the  Tsugaru  and  Soya  straits.  There  is 
not  much  information  about  the  volume  transport  through  the  Tatar  Strait  (Martin  and 
Kawase,  1998),  which  was  taken  as  0.05  Sv  in  this  study.  A  recent  estimate  of  the 
monthly  mean  volume  transport,  reported  by  Yi  (1966),  through  the  Korea/Tsushima 
Strait  with  the  annual  average  of  1.3  Sv,  a  maximum  of  2.2  Sv  in  October,  and  a 
minimum  of  0.3  Sv  in  February.  Bang  et  al.  (1996)  used  the  maximum  inflow  transport  of 
about  3.5  Sv  in  August  and  the  minimum  of  1.6  Sv  in  February.  In  contrast  to  Yi,  Kim 
and  Yoon  (1996)  used  the  mean  value  of  2.2  Sv  with  +/-  0.35  Sv  with  the  maximum  in 
mid-September  and  the  minimum  in  mid-March.  The  total  inflow  transport  through 
Korea/Tsushima  Straits  should  be  the  same  as  the  total  outflow  transport  through  the 
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Tsugaru  and  Soya  Straits.  We  assume  that  75%  (80%  in  Bang  et  al.,  1996)  of  the  total 
inflow  transport  should  flow  out  of  the  JES  through  the  Tsugaru  Strait,  and  25%  (20%  in 
Bang  et  al.,  1996)  through  the  Soya  Strait.  This  ratio  is  adopted  from  the  maximum 
volume  transport  through  the  Tsugaru  Strait  estimated  by  Toba  et  al.  (1982),  and  through 
the  Soya  Strait  estimated  by  Preller  and  Hogan  (1998). 


D.  MODE  SPLITTING 

For  computational  efficiency,  the  mode  splitting  technique  (Blumberg  and  Mellor, 
1987)  is  applied  with  a  barotropic  time  step  of  25  seconds,  based  on  the  Courant- 
Friederichs-Levy  computational  stability  (CFF)  condition  and  the  external  wave  speed; 
and  a  baroclinic  time  step  of  900  seconds,  based  on  the  CFF  condition  and  the  internal 
wave  speed. 


E.  TWO-STAGE  INITILIZATION 

1,  The  Pre-Simulation  Stage 

The  model  was  integrated  by  using  two-step  initialization.  During  the  first  step 
(restoring  run),  POM  is  integrated  for  two  years  from  zero  velocity  and  climatological 
temperature  and  salinity  fields  (Levitus,  1982)  with  the  climatological  monthly  mean 
surface  wind  stress  from  the  COADS  data  and  restoring-type  surface  thermohaline 
forcing  (ai  =  0,  012  =  1)  which  is  relaxed  to  surface  monthly  mean  values. 

2.  The  Simulation  Stage 

The  final  states  of  the  first  step  (restoring  run)  are  taken  as  initial  conditions  for 
the  second  step  (simulation  run).  During  the  simulation  run,  POM  is  integrated  again  for 
one  and  half  years  starting  from  Julian  Day  (JD)-l  to  JD-180  of  the  second  year  with 
monthly  mean  surface  wind  stress,  net  heat  flux,  and  fresh-water  flux  (ai  =  1,  012  =  0) 
from  the  COADS  data.  The  atmospheric  forcing  data  are  temporally  interpolated  into 
daily  data.  The  final  states  of  the  simulation  stage, 

Vo  =  Vjdiso,  To  =  Tjdi8o,  So  =  Sjdi80,  (6) 

are  taken  as  standard  initial  conditions  for  the  numerical  experiments. 
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IV.  THE  EXPERIMENTAL  DESIGN 


Uncertainty  in  the  initial/boundary  conditions,  and  atmospheric  forcing  leads  to 
the  uncertainty  in  model  output.  The  main  objective  of  this  study  is  to  investigate  the 
response  of  a  ocean  model  to  uncertain  input  data  using  POM  implemented  for  JES. 
Twelve  experiments  are  conducted  with  one  control  run  and  eleven  sensitivity  runs 
(Table  2). 


Experiment 

Property 

Description 

0 

Control  run 

See  section  A 

1 

Uncertain  velocity  initialization  processes 

See  section  B 

2 

3 

4 

5 

Uncertain  wind  stress 

See  section  C 

6 

7 

Uncertain  lateral  boundary  transport 

See  section  D 

8 

9 

Combination  of  uncertainty 

See  section  E 

10 

11 

Table  2  Summary  of  experimental  design. 
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A. 


THE  CONTROL  RUN 


The  control  run  is  to  integrate  POM-JES  from  the  standard  initial  conditions  (6) 
for  180  days  (to  JD-360)  with  the  lateral  transport  shown  in  Table  1  (unperturbed)  and 
the  daily  surface  wind  stress,  net  heat  flux,  and  fresh- water  flux  (ai  =  1,  ci2  =  0) 
interpolated  from  the  COADS  monthly  mean  data  (unperturbed).  The  simulated 
temperature  and  salinity  fields  and  circulation  pattern  are  consistent  with  observational 
studies  (Chu  et  al.  2001a). 


B.  UNCERTAIN  INITIAL  CONDITIONS 

As  mentioned  before,  initializing  the  velocity  field  with  the  diagnostic  mode 
(called  the  diagnostic  initialization)  contains  large  uncertainty  with  the  possibility  of 
generating  extremely  strong  thermohaline  source/sink  terms  (Chu  and  Lan,  2003).  Four 
experiments  are  designed  to  investigate  the  model  uncertainty  to  uncertain  initial  velocity 
conditions. 

Run-1  does  not  use  the  velocity  initialization.  The  POM-JES  prognostic  mode  is 
integrated  from 

Vo  =  0,  To  =  Tjdiso,  So  =  Sjdi80?  (7) 

with  the  same  atmospheric  and  lateral  boundary  forcing  as  Run-0  for  180  days.  Model 
difference  between  Run-0  and  Run- 1  is  the  uncertainty  caused  by  the  zero  initial  velocity 
fields. 

Run-2,  Run-3,  and  Run-4  are  designed  to  investigate  the  uncertainty  of  the 
diagnostic  initialization.  The  POM-JES  diagnostic  mode  is  integrated  from  (7)  (TJDigo, 
Sjdi80  unchanged)  for  30  days,  60  days,  and  90  days  to  get  \^g) ,  and  \^s) . 

The  POM-JES  prognostic  mode  is  integrated  with  the  same  atmospheric  and  lateral 
boundary  forcing  as  Run-0  for  180  days  from 

V0  =  V^\  To  =  Tjdiso,  So  =  Sjdi8o,  (8) 


in  Run-2;  from 
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V0=V™,  To=TJD180,  So  =  SJD180,  (9) 

in  Run-3;  and  from 

%=^90£)flg)’  To  =  Tjdi8o,  So  =  Sjdiso,  (10) 

in  Run-4  (Table  3). 


Experiment 

Initial  Conditions 

Wind  Forcing 

Lateral  Boundary 

Conditions 

1 

Vo=0, 

To=  TjD180, 

So  =  Sjdiso 

Same  as  Run-0 

Same  as  Run-0 

2 

AT  _  \r(Diag) 
y0  ~  y  30 D  ’ 

To  =TjD180, 

So  =  SjD180 

Same  as  Run-0 

Same  as  Run-0 

3 

AT  —  AT  (  Diag  ) 
y0  y  60  D  ’ 

To  =  TjD180, 

So  =  Sjd180 

Same  as  Run-0 

Same  as  Run-0 

4 

V  —  xr(Dias) 

y  0  y  90  D  5 

To  =  TjD180, 

So  =  Sjdiso 

Same  as  Run-0 

Same  as  Run-0 

Table  3  Experiments  for  uncertain  initial  conditions. 


C.  UNCERTAIN  WIND  FORCING 

Two  experiments  are  conducted  to  investigate  the  effect  of  wind  uncertainty. 
Everything  remains  the  same  as  Run-0  except  the  surface  winds  where  a  Gaussian-type 
random  variable  added  to  each  surface  grid  point  with  zero  mean  and  noise  intensity  of 
0.5  m  s"1  for  Run-5  and  1  m  s"1  for  Run-6,  respectively  (Table  4). 
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Experiment 

Initial  Conditions 

Wind  Forcing 

Lateral  Boundary 

Conditions 

5 

Same  as  Run-0 

Adding  Gaussian  random 

noise  with  zero  mean  and 

0.5  m/s  noise  intensity 

Same  as  Run-0 

6 

Same  as  Run-0 

Adding  Gaussian  random 

noise  with  zero  mean  and 

1.0  m/s  noise  intensity 

Same  as  Run-0 

Table  4  Experiments  for  uncertain  wind  forcing. 


D.  UNCERTAIN  LATERAL  TRANSPORT 

Two  experiments  are  conducted  to  investigate  the  effect  of  lateral  transport 
uncertainty.  Everything  keeps  the  same  as  Run-0  except  the  lateral  boundary  transport 
where  a  Gaussian-type  random  variable  is  added  to  each  lateral  boundary  grid  point  with 
the  zero  mean  and  noise  intensity  being  5%  and  10%  of  the  transport  (control  run)  for 
Run-7  and  Run-8,  respectively  (Table  5). 


Experiment 

Initial  Conditions 

Wind  Forcing 

Lateral  Boundary  Conditions 

7 

Same  as  Run-0 

Same  as  Run-0 

Adding  Gaussian  random  noise 

with  the  zero  mean  and  noise 

intensity  being  5%  of  the 

transport  (control  run) 

8 

Same  as  Run-0 

Same  as  Run-0 

Adding  Gaussian  random  noise 

with  the  zero  mean  and  noise 

intensity  being  10%  of  the 

transport  (control  run) 

Table  5  Experiments  for  uncertain  lateral  transport. 
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E.  COMBINED  UNCERTAINTY 

Three  experiments  are  conducted  to  investigate  the  effect  of  combined 
uncertainty.  Initial  conditions  remain  the  same  as  Run-2.  For  Run-9,  the  surface  wind  and 
the  lateral  boundary  transport  are  the  same  as  Run-6  and  Run-0,  respectively.  For  Run- 
10,  the  surface  wind  and  the  lateral  boundary  transport  are  the  same  as  Run-0  and  Run-8, 
respectively.  For  Run-11,  the  surface  wind  and  the  lateral  boundary  transport  are  the 
same  as  Run-6  and  Run-8,  respectively  (Table  6). 


Experiment 

Initial  conditions 

Wind  forcing 

Lateral  Boundary 

Conditions 

9 

Same  as  Run-2 

Same  as  Run-6 

Same  as  Run-0 

10 

Same  as  Run-2 

Same  as  Run-0 

Same  as  Run- 8 

11 

Same  as  Run-2 

Same  as  Run-6 

Same  as  Run- 8 

Table  6  Experiments  for  combined  uncertainty. 
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V.  STATISTICAL  ERROR  ANALYSIS 


A.  THE  MODEL  ERROR 

The  difference  between  the  variable  y/  o  f  control  run  and  each  experimental  run. 


A  y/{x,y,z,t)  =  y/c(x,y,z,t )  -  y/e{x,y,z,t) 


(ID 


represents  the  model  error.  Here  y/c  and  y/e  are  the  variables  (u  or  v  components  of  the 
horizontal  velocity  V)  from  control  run  and  each  experimental  run,  respectively.  This 
model  error  is  used  to  measure  the  strength  of  the  effect  on  each  experimental  run.  We 
may  take  the  vertically  averaging  value  and  probability  histogram  of  txy/  as  the  error 
distribution. 

B.  THE  ROOT  MEAN  SQUARE  ERROR 

The  Root  Mean  Square  Error  (RMSE)  of  the  horizontal  velocity  (V)  fields 
between  the  control  run  and  each  experimental  run  was  computed  as 


RMSE(z,t)  for  v 


1 


My  Mx 


\MyxMx  j=l  i=x 


,yj,z,tf+  A  y/v  (. Xi ,  yj,  z,  t)2  ] 


(12) 


where  Mx  and  M,  are  the  number  of  grid  points  along  the  east-west  direction  and  the 
north-south  direction  for  the  JES,  respectively. 

C.  THE  RELATIVE  ROOT  MEAN  SQUARE  ERROR 

The  Relative  Root  Mean  Square  Error  (RRMSE)  of  the  horizontal  velocity  (V) 
fields  between  the  control  run  and  each  experimental  run  were  also  computed  as 
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RRMSE(z,t ) 


(13) 


for  V 


My  IVlx 


y'=l  '=1 


My  Mi 


7=1  1=1 


+  Aif/Xxi,yj,z,t)2 


+  y/cv(xi,yj,z,tf 


where  Mx  and  M,  are  the  number  of  grid  points  along  the  east-west  direction  and  the 
north-south  direction  for  the  JES,  respectively. 
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VI.  MODEL  ERRORS  DUE  TO  INPUT  UNCERTAINTY 


A.  INITIAL  CONDITIONS 

1.  The  Model  Error 

Based  on  equation  (11)  derived  in  Chapter  V,  the  model  error  was  calculated.  For 
each  experimental  run  (Run  1  to  4),  the  u  and  v  components  of  the  horizontal  velocity  (V) 
fields  at  each  grid  point  were  subtracted  from  those  of  the  control  run  (Run  0).  These 
errors  were  used  to  measure  the  strength  of  the  effect  of  various  velocity  initialization 
processes  used  by  these  four  experimental  runs. 


a.  The  Horizontal  Distribution 

After  the  model  error  was  calculated  for  each  horizontal  grid  point,  the 
errors  for  all  the  22  sigma  levels  were  combined  vertically  for  each  field  on  the  5th  day 
(day-185)  and  the  180th  day  (day-360)  after  the  day  of  model  integration  (day-180).  Four 
plots  contained  the  horizontal  distribution  of  these  vertically  mean  model  errors  for  each 
specific  field  (a  and  v)  between  the  control  run  and  four  experimental  runs  made  on  the 
5th  day  and  the  180th  day.  The  values  of  max,  min,  mean  and  standard  deviation  were  also 
determined  for  each  plot.  These  plots  are  shown  in  Figure  4  to  7. 

As  seen  in  Figure  4,  for  model  errors  of  u  components  in  the  horizontal 
velocity  (V)  fields  on  the  5th  day,  numerous  small-  to  meso-scale  patterns  of  negative  and 
positive  model  errors  are  well-mixed  and  distributed  over  most  of  the  areas  of  JES.  Near 
the  East  Japan  Basin,  some  strong  negative  and  positive  patterns  are  present  in  Run  2  to 
4.  At  this  location,  both  the  maximum  negative  and  positive  model  errors  occur  in  Run  4 
with  a  value  of  -0.08  m/s  and  0.07  m/s,  respectively.  Some  isolated  strong  negative  and 
positive  model  errors  are  also  present  near  the  Ulleung  Basin  in  Run  1  with  a  value  of  - 
0.09  m/s  and  0.09  m/s,  respectively.  These  features  are  similar  for  the  model  errors  of  v 
components  in  the  horizontal  velocity  (V)  fields  on  the  5th  day  in  Figure  6.  Near  the  East 
Japan  Basin,  both  the  maximum  negative  and  positive  model  errors  occur  in  Run  4  with  a 

value  of  -0.13  m/s  and  0.13  m/s,  respectively.  Near  the  Ulleung  Basin,  some  isolated 
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strong  negative  and  positive  model  errors  are  also  present  in  Run  1  with  a  value  of  -0.08 
m/s  and  0.1  m/s,  respectively. 

As  seen  in  Figure  5,  model  errors  of  u  components  in  the  horizontal 
velocity  (V)  fields  on  the  180th  day  reveal  several  small-  to  meso-scale  patterns 
concentrated  near  the  Japan  Basin  on  each  plot.  Near  the  East  Japan  Basin,  some  strong 
negative  and  positive  patterns  are  present  for  all  the  experimental  runs.  At  this  location, 
the  maximum  negative  model  error  occurs  in  Run  2  with  a  value  of  -0.05  m/s,  which 
indicates  an  overestimation  in  model  prediction.  On  the  other  hand,  the  maximum 
positive  model  error  reached  0.03  °C  in  Run  2,  which  indicates  an  underestimation  in 
model  prediction.  These  features  are  similar  for  the  model  errors  of  v  components  in  the 
horizontal  velocity  (V)  fields  on  the  180th  day  in  Figure  7,  with  a  maximum  negative  and 
positive  value  of  -0.07  m/s  and  0.05  m/s  (Run  2),  respectively.  As  observed,  the  model 
errors  are  present  mostly  near  the  Japan  Basin.  This  could  be  associated  with  the  nearby 
Polar  Front  (PF)  where  two  current  systems  (NKCC  and  TWC)  converge.  In  general,  the 
model  errors  of  u  and  v  components  in  the  horizontal  velocity  (V)  fields  are  more 
significant  on  the  5th  day  than  those  on  the  180th  day.  In  addition,  not  much  difference 
exist  among  the  model  error  distributions  for  all  four  experimental  runs  on  both  days 
indicating  that  the  difference  among  the  model  errors  of  these  four  experimental  runs  is 
not  significant. 


b.  The  Histogram 

The  model  error,  calculated  by  equation  (11),  could  also  be  presented  in 
the  form  of  a  histogram  to  evaluate  the  model  error  distribution.  Four  histogram  plots 
contained  the  model  errors  in  each  specific  field  ( u  and  v)  between  the  control  run  and 
four  experimental  runs  made  on  the  5th  day  and  the  180th  day.  The  values  of  the  number 
of  the  total  samples,  samples  greater  than  zero,  samples  less  than  zero  and  standard 
deviation  were  also  determined  for  each  plot.  These  plots  are  shown  in  Figure  8  to  9. 

In  Figure  8,  the  histogram  of  model  errors  of  u  components  in  the 
horizontal  velocity  on  both  the  5th  day  and  180th  day  reveals  a  Gaussian-type  distribution 
on  each  plot.  The  values  of  standard  deviation  (STD)  are  decreasing  with  integration 
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time.  For  the  5th  day,  the  maximum  and  minimum  values  of  standard  deviation  occur  in 
Run  1  and  Run  2  with  a  value  of  0.023  m/s  and  0.021  m/s,  respectively.  For  the  180th  day, 
the  maximum  and  minimum  values  of  standard  deviation  occur  in  Run  1  and  Run  3  with 
a  much  smaller  value  of  0.008  m/s  and  0.007  m/s,  respectively.  These  features  are  similar 
for  the  histogram  of  model  errors  of  v  components  in  the  horizontal  velocity  on  both  the 
5th  day  and  180th  day  in  Figure  9.  For  the  5th  day,  the  maximum  and  minimum  values  of 
standard  deviation  occur  in  Run  4  and  Run  2  with  a  value  of  0.025  m/s  and  0.023  m/s, 
respectively.  For  the  180th  day,  the  maximum  and  minimum  values  of  standard  deviation 
are  occurred  in  Run  1  and  Run  4  with  a  value  of  0.009  m/s  and  0.008  m/s,  respectively.  In 
general,  obvious  differences  exist  between  these  two  days  (the  model  errors  are 
decreasing  with  time),  but  the  differences  among  these  four  experimental  runs  are  still  not 
significant. 
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Latitude  Latitude 


50°N 


40°N 


45°N 


(a) 


Max  :  0.085936 
Min  :  -0.090132 
Mean  :  0.0011568 
Std  :  0.017027 


(b) 


Max  :  0.048645 
Min  :  -0.057586 
Mean  :  0.001112 
Std  :  0.013917 


Figure  4.  Vertically  averaged  model  error  of  (a)  Run  0  -  Run  1,  (b)  Run  0  -  Run  2, 
(c)  Run  0  -  Run  3  and  (d)  Run  0  -  Run  4  for  u  component  of  velocity  field 
on  the  5th  day  after  the  model  integration. 
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Latitude  Latitude 


Figure  5.  Vertically  averaged  model  error  of  (a)  Run  0  -  Run  1,  (b)  Run  0  -  Run  2, 
(c)  Run  0  -  Run  3  and  (d)  Run  0  -  Run  4  for  u  component  of  velocity  field 
on  the  180th  day  after  the  model  integration. 
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Latitude  Latitude 


50°N 


40°N 
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Figure  6.  Vertically  averaged  model  error  of  (a)  Run  0  -  Run  1,  (b)  Run  0  -  Run  2, 
(c)  Run  0  -  Run  3  and  (d)  Run  0  -  Run  4  for  v  component  of  velocity  field 
on  the  5th  day  after  the  model  integration. 
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Figure  7.  Vertically  averaged  model  error  of  (a)  Run  0  -  Run  1,  (b)  Run  0  -  Run  2, 
(c)  Run  0  -  Run  3  and  (d)  Run  0  -  Run  4  for  v  component  of  velocity  field 
on  the  180th  day  after  the  model  integration. 
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Figure  8.  Histogram  of  (a)  Run  0  -  Run  1,  (b)  Run  0  -  Run  2,  (c)  Run  0  -  Run  3  and 
(d)  Run  0  -  Run  4  for  u  component  of  velocity  field  on  both  the  5th  and 
180th  day  after  the  model  integration. 
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Figure  9.  Histogram  of  (a)  Run  0  -  Run  l,  (b)  Run  0  -  Run  2,  (c)  Run  0  -  Run  3  and 
(d)  Run  0  -  Run  4  for  v  component  of  velocity  field  on  both  the  5th  and 
180th  day  after  the  model  integration. 
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2. 


The  Root  Mean  Square  Error  (RMSE) 


a.  The  Vertical  Variation 

The  Root  Mean  Square  Error  (RMSE),  calculated  by  equation  (12),  is 
commonly  used  to  evaluate  the  model  performance.  For  the  vertical  profile  of  RMSE, 
two  plots  contained  the  errors  caused  by  these  four  experimental  runs  for  the  horizontal 
velocity  (V)  fields  made  on  the  both  the  5th  day  (day- 185)  and  the  180th  day  (day-360) 
after  the  day  of  model  integration  (day- 180).  These  plots  are  shown  in  Figure  10. 

As  seen  in  this  figure,  RMSE  on  the  5th  day  reveals  higher  value  above 
sigma  level  8,  which  indicates  some  stronger  effects  on  model  prediction  occurring  at 
those  levels.  The  maximum  RMSE  occurs  in  Run  1  with  a  value  of  0.056  m/s  at  level  1. 
These  features  are  similar  for  the  RMSE  on  the  180th  day.  The  maximum  RMSE  occurs 
in  Run  2  with  a  value  of  0.02  m/s  at  level  1.  In  general,  from  these  plots,  the  curves 
representing  these  four  experimental  runs  are  very  close  to  each  other,  indicating  no 
obvious  difference  among  the  effects  of  these  four  experimental  runs.  In  addition,  the 
RMSE  of  horizontal  velocity  (V)  changed  dramatically  from  the  5th  day  (0.05  m/s  near 
surface  level  and  0.03  m/s  near  bottom  level)  to  the  180th  day  (0.02  m/s  near  surface  level 
and  0.01  m/s  near  bottom  level). 

b.  The  Temporal  Evolution 

The  temporal  evolution  of  vertically  averaged  RMSE  for  the  horizontal 
velocity  (V)  fields  from  the  5th  day  to  180th  day  is  shown  in  Figure  11.  As  seen  in  this 
figure,  no  significant  difference  shows  among  these  four  curves  of  RMSE  for  the 
horizontal  velocity  (V)  fields.  The  RMSE  of  horizontal  velocity  rapidly  decreases  with 
time  in  the  first  20  days  from  a  peak  value  of  0.032  m/s  to  0.02  m/s.  It  then  slowly 
decreases  with  time  to  0.01  m/s  on  the  180th  day.  In  general,  as  in  the  previous  results,  no 
obvious  difference  exists  among  the  effects  of  these  four  experimental  runs. 
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Figure  10.  RMSE  of  V  in  (a)  Run  0  -  Run  1,  (b)  Run  0  -  Run  2,  (c)  Run  0  -  Run  3 
and  (d)  Run  0  -  Run  4  on  the  5th  day  and  the  180th  day  after  the  model 
integration. 


Figure  11.  Temporal  evolution  of  vertically  averaged  RMSE  of  V  in  a)  Run  0  -  Run 
1 ,  (b)  Run  0  -  Run  2,  (c)  Run  0  -  Run  3  and  (d)  Run  0  -  Run  4. 
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3. 


The  Relative  Root  Mean  Square  Error  (RRMSE) 


a.  The  Vertical  Variation 

The  Relative  Root  Mean  Square  Error  (RRMSE),  calculated  by  equation 
(13),  is  also  commonly  used  for  evaluating  model  performance.  For  the  vertical  profile  of 
RRMSE,  two  plots  contained  the  errors  caused  by  these  four  experimental  runs  for  the 
horizontal  velocity  (V)  fields  made  on  both  the  5th  day  (day- 185)  and  the  180th  day  (day- 
360)  after  the  day  of  model  integration  (day-180).  These  plots  are  shown  in  Figure  12. 

As  seen  in  this  figure,  the  RRMSE  on  the  5  th  day  reveals  a  higher  value 
below  sigma  level  18,  which  indicates  some  stronger  effects  on  model  prediction 
occurring  at  those  levels.  The  maximum  RRMSE  occurs  in  Run  1  with  a  value  of  0.75  at 
level  22.  On  the  other  hand,  the  RRMSE  on  the  180th  day  reveals  a  higher  value  above 
sigma  level  8,  which  indicates  some  stronger  effects  on  model  prediction  occurring  at 
those  levels.  The  maximum  RRMSE  occurs  in  Run  2  with  a  value  of  0.26  at  level  4.  In 
general,  from  these  plots,  the  curves  representing  these  four  experimental  runs  are  very 
close  to  each  other,  indicating  no  obvious  difference  among  the  effects  of  these  four 
experimental  runs.  Similar  to  previous  plots  of  the  RMSE,  the  RRMSE  of  horizontal 
velocity  (V)  changed  dramatically  from  the  5th  day  (0.7  near  surface  level  and  0.5  near 
bottom  level)  to  the  180th  day  (0.25  near  surface  level  and  0.2  near  bottom  level). 

b.  The  Temporal  Evolution 

The  temporal  evolution  of  vertically  averaged  RRMSE  for  the  horizontal 
velocity  (V)  fields  from  the  5th  day  to  the  180th  day  is  shown  in  Figure  13.  As  seen  in  this 
figure,  no  significant  difference  exists  among  these  four  curves  of  the  RRMSE  for  the 
horizontal  velocity  (V)  fields.  The  RRMSE  of  horizontal  velocity  rapidly  decreases  with 
time  in  the  first  20  days  from  a  peak  value  of  0.5  to  0.3  and  then  slowly  decreases  with 
time  to  0.2  on  the  180th  day.  In  general,  no  obvious  difference  among  the  effects  of  these 
four  experimental  runs  exists.  Notably,  the  RRMSE  of  horizontal  velocity  was  between 
0.3  and  0.5.  Therefore,  the  effects  of  the  velocity  initialization  processes  to  the  model 
horizontal  velocity  prediction  are  quite  significant. 
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Figure  12.  RRMSE  of  V  in  (a)  Run  0  -  Run  l,  (b)  Run  0  -  Run  2,  (c)  Run  0  -  Run  3 
and  (d)  Run  0  -  Run  4  on  the  5th  day  and  the  180th  day  after  the  model 
integration. 


Figure  13.  Temporal  evolution  of  vertically  averaged  RRMSE  of  V  in  (a)  Run  0  - 
Run  1 ,  (b)  Run  0  -  Run  2,  (c)  Run  0  -  Run  3  and  (d)  Run  0  -  Run  4. 
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B.  WIND  FORCING 

1.  The  Model  Error 

a.  The  Horizontal  Distribution 

As  mentioned  in  section  A,  after  the  model  error  was  calculated,  for  each 
horizontal  grid  point,  the  errors  for  all  the  22  sigma  levels  were  combined  vertically  for 
each  field  on  the  5th  day  (day- 185)  and  the  180th  day  (day-360)  after  the  day  of  model 
integration  (day- 180).  Four  plots  contained  the  horizontal  distribution  of  these  vertically 
mean  model  errors  for  each  specific  field  (u  and  v)  between  the  control  run  and  two 
experimental  runs  (Run  5  and  Run  6)  made  on  both  the  5th  day  and  180th  day.  The  values 
of  max,  min,  mean  and  standard  deviation  were  also  determined  for  each  plot.  These 
plots  are  shown  in  Figure  14  to  15. 

As  seen  in  these  figures,  the  model  errors  of  both  u  and  v  on  the  5th  and 
180th  days  reveal  that  numerous  small-  to  meso-scale  patterns  of  negative  and  positive 
model  errors  are  distributed  near  the  Japan  Basin  and  some  of  the  coastal  areas  of  JES. 
As  expected,  for  the  model  errors  on  the  5th  day,  relatively  strong  patterns  are  present  in 
Run  6  with  relatively  high  values  of  the  maximum,  minimum,  mean  and  standard 
deviation.  A  significant  difference  in  these  statistic  values  exists  between  the  two 
experimental  runs.  These  features  are  similar  for  the  model  errors  on  the  180th  day.  These 
relatively  strong  patterns  are  also  present  in  Run  6  with  the  higher  statistic  values.  This 
indicates  the  model  error  caused  by  Run  6  had  a  larger  effect  on  the  model  performance 
than  that  caused  by  the  Run  5.  As  observed,  the  model  errors  are  present  mostly  near  the 
Japan  Basin.  This  could  be  associated  with  the  nearby  Polar  Front  (PF)  where  two  current 
systems  (NKCC  and  TWC)  converge.  In  general,  from  these  horizontal  distribution  plots, 
obvious  differences  exist  between  Run  5  and  Run  6  (more  model  error  was  caused  by 
Run  6).  In  addition,  the  model  errors  increase  slowly  from  the  5th  to  the  180th  day. 
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b. 


The  Histogram 


The  model  error,  calculated  by  equation  (11),  is  also  present  in  the  form  of 
a  histogram  in  evaluating  the  model  error  distribution.  Two  histogram  plots  contained  the 
model  errors  on  each  specific  field  ( u  and  v)  between  the  control  run  and  two 
experimental  runs  made  on  the  5th  day  and  the  180th  day.  The  values  of  the  number  of  the 
total  sample,  samples  greater  than  zero,  samples  less  than  zero  and  standard  deviation 
were  also  determined  for  each  plot.  These  plots  are  shown  in  Figure  16  to  17. 

In  Figure  16,  the  histogram  of  model  errors  of  u  components  of  the 
horizontal  velocity  on  both  the  5th  day  and  180th  day  reveals  a  Gaussian-type  distribution 
on  each  plot.  Similar  to  previous  plots,  the  values  of  standard  deviation  (STD)  are  higher 
in  Run  6  than  those  in  Run  5  and  they  increase  with  integration  time  on  each 
experimental  run.  For  the  5  th  day,  the  maximum  and  minimum  values  of  standard 
deviation  occur  in  Run  6  and  Run  5  with  a  value  of  0.015  m/s  and  0.009  m/s, 
respectively.  For  the  180th  day,  the  maximum  and  minimum  values  of  standard  deviation 
also  occur  in  Run  6  and  Run  5  with  a  value  of  0.016  m/s  and  0.01  m/s,  respectively. 
These  features  are  similar  for  the  histogram  of  model  errors  of  v  components  of  the 
horizontal  velocity  on  both  the  5th  day  and  180th  day  in  Figure  17.  Again,  for  both  the  5th 
and  180th  day,  the  maximum  and  minimum  values  of  standard  deviation  occur  in  Run  6 
and  Run  5  with  the  maximum  and  minimum  value  of  0.017  m/s  and  0.001  m/s  on  the  5th 
day,  and  0.018  m/s  and  0.001  m/s  on  the  180th  day,  respectively.  Generally  speaking,  the 
differences  between  Run  5  and  Run  6  (more  model  error  was  caused  by  Run  6)  are  quite 
significant  while  these  model  errors  increase  slightly  from  the  5th  day  to  the  180th  day. 
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Figure  14.  Vertically  averaged  model  error  of  (a)  Run  0  -  Run  5  and  (b)  Run  0  -  Run 
6  for  u  component  of  velocity  field  on  both  the  5th  day  and  180th  day  after 
the  model  integration. 
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Figure  15.  Vertically  averaged  model  error  of  (a)  Run  0  -  Run  5  and  (b)  Run  0  -  Run 
6  for  v  component  of  velocity  field  on  both  the  5th  day  and  180th  day  after 
the  model  integration. 
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Figure  16. 


Histogram  of  (a)  Run  0  -  Run  5  and  (b)  Run  0  -  Run  6  for  u  component  of 
velocity  field  on  both  the  5th  and  180th  day  after  the  model  integration. 
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Figure  17. 


Histogram  of  (a)  Run  0  -  Run  5  and  (b)  Run  0  -  Run  6  for  v  component  of 
velocity  field  on  both  the  5th  and  180th  day  after  the  model  integration. 
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2. 


The  Root  Mean  Square  Error  (RMSE) 


a.  The  Vertical  Variation 

The  Root  Mean  Square  Error  calculated  by  equation  (12)  evaluates  the 
model  performance.  For  the  vertical  profile  of  RMSE,  two  plots  contained  the  errors 
caused  by  these  two  experimental  runs  for  the  horizontal  velocity  (V)  fields  made  on  both 
the  5th  day  (day- 185)  and  the  180th  day  (day-360)  after  the  day  of  model  integration  (day- 
180).  These  plots  are  shown  in  Figure  18. 

As  seen  in  this  figure,  the  RMSE  on  the  5th  day  reveals  a  higher  value 
occurring  above  sigma  level  8,  which  indicates  that  some  stronger  effects  on  model 
prediction  occurring  at  those  levels.  The  maximum  RMSE  occurs  in  Run  6  at  level  1  with 
a  value  of  0.068  m/s.  These  features  are  similar  for  the  RMSE  on  the  180th  day.  As 
expected,  the  maximum  RMSE  of  horizontal  velocity  occurs  in  Run  6  at  level  1  with  a 
value  of  0.07  m/s.  In  general,  from  these  plots,  obvious  differences  exist  between  the 
RMSE  in  Run  5  and  Run  6  (more  RMSE  was  caused  by  Run  6)  increasing  slowly  from 
the  5th  day  to  the  180th  day.  In  addition,  these  values  of  RMSE  decrease  with  depth 
(rapidly  above  sigma  level  8). 

b.  The  Temporal  Evolution 

The  temporal  evolution  of  vertically  averaged  RMSE  for  horizontal 
velocity  (V)  fields  from  the  5th  day  to  180th  day  is  shown  in  Figure  19.  As  the  figure 
demonstrates,  significant  differences  among  these  two  curves  of  RMSE  are  evident  for 
the  horizontal  velocity  (V)  fields  (more  RMSE  was  caused  by  Run  6).  The  RMSE 
increases  with  time  in  the  first  45  days  and  reaches  a  peak  value  of  0.018  m/s  and  0.01 1 
m/s  in  Run  6  and  Run  5,  respectively.  It  then  rapidly  decreases  from  the  45th  day  to  the 
60th  day.  The  changes  are  little  between  the  60th  day  and  the  130th  day,  but  they  increase 
again  from  the  130th  day  to  the  180th  day.  For  Run  6  and  Run  5,  the  increased  value  from 
the  5th  day  to  the  180th  day  is  0.003  m/s  and  0.002  m/s,  respectively.  In  general,  just  like 
the  other  previous  results,  obvious  differences  occur  between  the  RMSE  in  Run  5  and 
Run  6  (more  RMSE  was  caused  by  Run  6)  increasing  from  the  5th  day  to  the  180th  day. 
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Figure  18.  RMSE  of  V  in  (a)  Run  0  -  Run  5  and  (b)  Run  0  -  Run  6  on  the  5th  day  and 
the  180th  day  after  the  model  integration. 


Figure  19.  Temporal  evolution  of  vertically  averaged  RMSE  of  V  in  (a)  Run  0  -  Run 
5  and  (b)  Run  0  -  Run  6. 
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3. 


The  Relative  Root  Mean  Square  Error  (RRMSE) 


a.  The  Vertical  Variation 

The  Relative  Root  Mean  Square  Error,  calculated  by  equation  (13),  is  also 
commonly  used  when  evaluating  the  model  performance.  For  the  vertical  profile  of 
RRMSE,  two  plots  contained  the  errors  caused  by  these  two  experimental  runs  for  the 
horizontal  velocity  (V)  fields  made  on  both  the  5th  day  (day- 185)  and  the  180th  day  (day- 
360)  after  the  day  of  model  integration  (day-180).  These  plots  are  shown  in  Figure  20. 

As  seen  in  this  figure,  the  RRMSE  on  the  5  th  day  reveals  a  higher  value 
above  sigma  level  8,  which  indicates  some  stronger  effects  on  model  prediction  occurring 
at  those  levels.  The  maximum  RRMSE  occurs  in  Run  6  with  a  value  of  0.58  at  level  1  for 
horizontal  velocity.  The  RRMSE  on  the  180th  day  also  reveals  a  higher  value  above 
sigma  level  8,  which  indicates  some  stronger  effects  on  model  prediction  occurring  at 
those  levels.  As  expected,  the  maximum  RRMSE  occurs  in  Run  6  with  a  value  of  0.76  at 
level  1  for  horizontal  velocity.  In  general,  from  these  plots,  the  RRMSE  caused  by  Run  6 
is  larger  than  that  caused  by  Run  5  on  both  days.  Generally,  the  RRMSE  decreases  with 
depth  (rapidly  above  sigma  level  8).  As  observed,  the  effect  of  the  wind  forcing  noise  on 
the  model  velocity  field  prediction  is  significant,  especially  above  sigma  level  8  (around 
60  %  and  35  %  on  the  5th  day  and  80  %  and  50  %  on  the  180th  day  at  the  surface  level  in 
Run  6  and  Run  5,  respectively). 

b.  The  Temporal  Evolution 

The  temporal  evolution  of  vertically  averaged  RRMSE  for  horizontal 
velocity  (V)  fields  from  the  5th  day  to  the  180th  day  is  shown  in  Figure  21.  As  seen  in  this 
figure,  significant  differences  exist  between  these  two  curves  of  RRMSE  for  the 
horizontal  velocity  (V)  fields  (more  RRMSE  was  caused  by  Run  6).  The  RRMSE 
increases  with  time  slowly  in  the  first  45  days  and  then  decreases  rapidly  from  the  45th 
day  to  the  60th  day.  These  changes  are  little  between  the  60th  day  and  the  120th  day  but 
then  increase  again  from  the  130th  day  to  the  180th  day.  The  maximum  and  minimum 
RRMSE  of  horizontal  velocity  is  0.28  and  0.12  in  Run  6  and  0.19  and  0.08  in  Run  5, 
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respectively.  As  in  previous  results,  obvious  differences  take  place  between  the  RRMSE 
in  Run  5  and  Run  6  (more  RRMSE  was  caused  by  Run  6)  increasing  from  the  5th  day  to 
the  180th  day. 
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Figure  20.  RRMSE  of  V  in  (a)  Run  0  -  Run  5  and  (b)  Run  0  -  Run  6  on  the  5th  day 
and  the  180th  day  after  the  model  integration. 


Figure  21.  Temporal  evolution  of  vertically  averaged  RRMSE  of  V  in  (a)  Run  0  - 
Run  5  and  (b)  Run  0  -  Run  6. 
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c. 


OPEN  BOUNDARY  CONDITIONS 


1.  The  Model  Error 

a.  The  Horizontal  Distribution 

As  mentioned  in  section  A,  after  the  model  error  was  calculated,  for  each 
horizontal  grid  point,  the  errors  for  all  the  22  sigma  levels  were  combined  vertically  for 
each  field  on  the  5th  day  (day- 185)  and  the  180th  day  (day-360)  after  the  day  of  model 
integration  (day- 180).  Four  plots  contained  the  horizontal  distribution  of  these  vertically 
mean  model  errors  in  each  specific  field  (u  and  v)  between  the  control  run  and  two 
experimental  runs  (Run  7  and  Run  8)  made  on  the  5th  day  and  the  180th  day.  The  values 
of  max,  min,  mean  and  standard  deviation  were  also  determined  for  each  plot.  These 
plots  are  shown  in  Figure  22  to  23. 

As  seen  in  these  figures,  model  errors  of  both  u  and  v  on  the  5th  and  the 
180th  day  reveal  that  only  a  few  small-scale  patterns  occur  near  the  Tatar  Strait  and  the 
Soya  Strait  on  the  5th  day;  however,  some  small-  to  meso-scale  patterns  of  model  errors 
are  present  near  the  Tatar  Strait  and  the  Japan  Basin  on  the  180th  day.  For  the  model 
errors  on  the  5  th  day,  relatively  strong  patterns  are  present  in  Run  8  with  relatively  high 
values  of  the  maximum,  minimum,  mean  and  standard  deviation.  These  features  are 
similar  for  the  model  errors  on  the  180th  day  with  increased  statistic  values.  These 
relatively  strong  patterns  are  also  present  in  Run  8  with  the  higher  statistic  values.  This 
indicates  that  the  model  error  caused  by  Run  8  had  a  larger  effect  on  the  model 
performance  than  that  caused  by  Run  7.  In  general,  from  these  horizontal  distribution 
plots,  obvious  differences  between  Run  7  and  Run  8  exist  (more  model  error  was  caused 
by  Run  8)  with  the  model  errors  also  increasing  from  the  5th  to  the  180th  day. 

b.  The  Histogram 

The  model  error,  calculated  by  equation  (11),  is  also  present  in  the  form  of 

a  histogram  when  evaluating  the  model  error  distribution.  Two  histogram  plots  contained 

the  model  errors  on  each  specific  field  (it  and  v)  between  the  control  run  and  two 
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experimental  runs  made  on  the  5th  day  and  the  180th  day.  The  values  of  the  number  of  the 
total  sample,  samples  greater  than  zero,  samples  less  than  zero  and  standard  deviation 
were  also  determined  for  each  plot.  These  plots  are  shown  in  Figure  24  to  25. 

In  Figure  24,  the  histogram  of  model  errors  of  u  components  of  the 
horizontal  velocity  on  both  the  5th  day  and  180th  day  reveals  a  Gaussian-type  distribution 
on  each  plot.  Similar  to  the  previous  plots,  the  values  of  standard  deviation  (STD)  are 
higher  in  Run  8  than  those  in  Run  7  increasing  with  integration  time  on  each 
experimental  run.  For  the  5  th  day,  the  maximum  and  minimum  values  of  standard 
deviation  occur  in  Run  8  and  Run  7  with  a  value  of  0.006  m/s  and  0.003  m/s, 
respectively.  For  the  180th  day,  the  maximum  and  minimum  values  of  standard  deviation 
also  occur  in  Run  8  and  Run  7  with  a  value  of  0.008  m/s  and  0.005  m/s,  respectively. 
These  features  are  similar  for  the  histogram  of  model  errors  of  v  components  of  the 
horizontal  velocity  on  both  the  5th  day  and  the  180th  day  in  Figure  25.  Once  again,  for 
both  the  5th  and  the  180th  day,  the  maximum  and  minimum  values  of  standard  deviation 
occur  in  Run  8  and  Run  7  with  a  value  of  0.007  m/s  and  0.004  m/s  on  the  5th  day,  and 
0.009  m/s  and  0.006  m/s  on  the  180th  day,  respectively.  Generally  speaking,  these  model 
errors  increase  from  the  5th  day  to  the  180th  day  and  the  differences  between  Run  7  and 
Run  8  (more  model  error  was  caused  by  Run  8)  are  significant. 
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Figure  22.  Vertically  averaged  model  error  of  (a)  Run  0  -  Run  7  and  (b)  Run  0  -  Run 
8  for  u  component  of  velocity  field  on  both  the  5th  day  and  180th  day  after 
the  model  integration. 


50 


Latitude  Latitude 


180th  Day 


5th  Day 


Figure  23.  Vertically  averaged  model  error  of  (a)  Run  0  -  Run  7  and  (b)  Run  0  -  Run 
8  for  v  component  of  velocity  field  on  both  the  5th  day  and  180th  day  after 
the  model  integration. 
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Figure  24.  Histogram  of  (a)  Run  0  -  Run  7  and  (b)  Run  0  -  Run  8  for  u  component  of 
velocity  field  on  both  the  5th  and  180th  day  after  the  model  integration. 
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Figure  25.  Histogram  of  (a)  Run  0  -  Run  7  and  (b)  Run  0  -  Run  8  for  v  component  of 
velocity  field  on  both  the  5th  and  180th  day  after  the  model  integration. 
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2. 


The  Root  Mean  Square  Error  (RMSE) 


a.  The  Vertical  Variation 

The  Root  Mean  Square  Error,  calculated  by  equation  (12),  is  used  for 
evaluating  the  model  performance.  For  the  vertical  profile  of  RMSE,  two  plots  contained 
the  errors  caused  by  these  two  experimental  runs  for  the  horizontal  velocity  (V)  fields 
made  on  both  the  5th  day  (day- 185)  and  the  180th  day  (day-360)  after  the  day  of  model 
integration  (day- 180).  These  plots  are  shown  in  Figure  26. 

As  seen  in  this  figure,  the  RMSE  on  the  5th  day  reveals  that  the  higher 
value  of  RMSE  is  caused  by  Run  8  on  both  days.  The  maximum  RMSE  occurs  in  Run  8 
with  a  value  of  0.001  m/s  at  level  21.  These  features  are  similar  for  the  RMSE  on  the 
180th  day.  However,  the  relatively  higher  value  shifts  upward  to  above  sigma  level  8, 
which  indicates  that  some  stronger  effects  on  model  prediction  occurs  at  those  levels.  As 
expected,  the  maximum  RMSE  occurs  in  Run  8  at  level  1  with  a  value  of  0.015  m/s.  In 
general,  from  these  plots,  obvious  differences  exist  between  the  RMSE  in  Run  7  and  Run 
8  (more  RMSE  are  caused  by  Run  8)  increasing  from  the  5th  day  to  the  180th  day 
(obviously  above  the  sigma  level  8). 

b.  The  Temporal  Evolution 

The  temporal  evolution  of  vertically  averaged  RMSE  for  horizontal 
velocity  (V)  fields  from  the  5th  day  to  the  180th  day  is  shown  in  Figure  27.  As  seen  in  this 
figure,  a  significant  difference  exists  between  these  two  curves  of  RMSE  for  the 
horizontal  velocity  (V)  fields  (more  RMSE  is  caused  by  Run  8).  The  RMSE  of  horizontal 
velocity  oscillates  and  generally  increases  with  time  from  the  5th  day  to  the  90th  day  with 
a  peak  value  of  0.022  m/s  and  0.012  m/s  in  Run  8  and  Run  7,  respectively.  It  then 
oscillates  and  generally  decreases  with  time  from  the  90th  day  to  the  180th  day.  For  both 
Run  8  and  Run  7,  the  increased  value  from  the  5th  day  to  the  180th  day  is  0.003  m/s.  As  in 
previous  results,  obvious  differences  exist  between  the  RMSE  in  Run  7  and  Run  8  (more 
RMSE  is  caused  by  Run  8)  also  increasing  from  the  5th  day  to  the  180th  day. 


53 


RMSE  of  V  (m/s) 


180th  Day 


5th  Day 


O 

CO 


2  - 

4  - 

6 

8 

10 

12  h 

14 

16 

18 

20 

22 


AC) 

O 

2i  O 

2^  o 
^  o 
2i  o 
2i  o 
2S.  o 

2S.  O 
2^  O 
•2 i  o 
A  o 

A() 
2^  O 
AC) 
2i  O 
2VC> 
2i  O 
•  2i  cb  ■ 

2i  6 


0.02 

RMSE  of  V  (m/s) 


0.04 


Figure  26.  RMSE  of  V  in  (a)  Run  0  -  Run  7  and  (b)  Run  0  -  Run  8  on  the  5th  day  and 
the  180th  day  after  the  model  integration. 


Figure  27.  Temporal  evolution  of  vertically  averaged  RMSE  of  V  in  (a)  Run  0  -  Run 
7  and  (b)  Run  0  -  Run  8. 
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3. 


The  Relative  Root  Mean  Square  Error  (RRMSE) 


a.  The  Vertical  Variation 

The  Relative  Root  Mean  Square  Error,  calculated  by  equation  (13),  is  also 
commonly  used  for  evaluating  the  model  performance.  For  the  vertical  profile  of 
RRMSE,  two  plots  contained  the  errors  caused  by  these  two  experimental  runs  for  the 
horizontal  velocity  (V)  fields  made  on  both  the  5th  day  (day- 185)  and  the  180th  day  (day- 
360)  after  the  day  of  model  integration  (day-180).  These  plots  are  shown  in  Figure  28. 

As  seen  in  this  figure,  the  RRMSE  on  the  5th  day  reveals  that  a  relatively 
higher  value  occurs  below  sigma  level  8,  which  indicates  some  stronger  effects  on  model 
prediction  occurring  at  those  levels.  The  higher  value  of  RRMSE  is  caused  by  Run  8  on 
both  days  (the  deeper  the  level,  the  higher  the  value  of  difference).  The  maximum 
RRMSE  occurs  in  Run  8  with  a  value  of  0.23  at  level  22.  These  features  are  similar  for 
the  RRMSE  on  the  180th  day.  The  relatively  higher  RRMSE  occurs  in  Run  8.  As 
expected,  the  maximum  RRMSE  occurs  in  Run  8  with  a  value  of  0  0.28  at  level  22.  In 
general,  in  these  plots,  the  RRMSE  caused  by  Run  8  is  larger  than  that  by  Run  7  on  both 
days  and  increases  with  depth.  As  observed,  the  effects  of  the  uncertainty  on  lateral 
boundary  transport  to  the  model  velocity  field  prediction  (8  %  at  level  1  and  25  %  at  level 
22  on  the  5th  day  and  18  %  at  level  1  and  28  %  at  level  22  on  the  180th  day  in  Run  8)  are 
also  significant. 

b.  The  Temporal  Evolution 

The  temporal  evolution  of  vertically  averaged  RRMSE  for  horizontal 
velocity  (V)  fields  from  the  5th  day  to  180th  day  is  shown  in  Figure  29.  As  seen  from  this 
figure,  a  significant  difference  exists  between  these  two  curves  of  RRMSE  for  the 
horizontal  velocity  (V)  fields  (more  RRMSE  is  caused  by  Run  8).  The  RRMSE  of 
horizontal  velocity  oscillates  and  generally  increases  with  time  from  the  5th  day  to  the 
150th  day  with  a  peak  value  of  0.32  and  0.2  in  Run  8  and  Run  7,  respectively.  It  then 
oscillates  and  generally  decreases  with  time  from  the  150th  day  to  the  180th  day.  For  both 
Run  8  and  Run  7,  the  increased  value  from  the  5th  day  to  the  180th  day  is  0.05.  As  in 
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previous  results,  obvious  differences  exist  between  the  RRMSE  in  Run  7  and  Run  8 
(more  RRMSE  is  caused  by  Run  8),  also  increasing  from  the  5th  day  to  the  180th  day. 
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Figure  28.  RRMSE  of  V  in  (a)  Run  0  -  Run  7  and  (b)  Run  0  -  Run  8  on  the  5th  day 
and  the  180th  day  after  the  model  integration. 


Figure  29.  Temporal  evolution  of  vertically  averaged  RRMSE  of  V  in  (a)  Run  0  - 
Run  7  and  (b)  Run  0  -  Run  8. 
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D.  COMBINED  UNCERTAINTY 


1.  The  Model  Error 

a.  The  Horizontal  Distribution 

Based  on  (11)  derived  in  Chapter  VI,  the  model  error  was  calculated. 
These  errors  are  used  to  measure  the  strength  of  the  effect  of  various  combined 
uncertainties  used  by  these  three  experimental  runs.  After  calculating  the  model  error,  for 
each  horizontal  grid  point,  the  errors  for  all  the  22  sigma  levels  were  combined  vertically 
for  each  field  on  the  5th  day  (day- 185)  and  the  180th  day  (day-360)  after  the  day  of  model 
integration  (day- 180).  Six  plots  contained  the  horizontal  distribution  of  these  vertically 
mean  model  errors  on  each  specific  field  (u  and  v)  between  the  control  run  and  three 
experimental  runs  (Run  9,  10  and  11)  made  on  the  5th  day  and  the  180th  day.  The  values 
of  max,  min,  mean  and  standard  deviation  were  also  determined  for  each  plot.  These 
plots  are  shown  in  Figure  30  to  3 1 . 

As  seen  in  these  figures,  the  model  errors  of  both  u  and  v  on  the  5th  and 
the  180th  days  reveal  that  numerous  small-  to  meso-scale  patterns  of  negative  and  positive 
model  errors  are  distributed  near  the  Japan  Basin  and  the  Tatar  Strait  (especially  on  the 
180th  day).  For  the  model  errors  on  the  5th  day,  as  expected,  relatively  strong  patterns  are 
present  in  Run  1 1  with  relatively  high  statistic  values  of  the  maximum,  minimum,  mean 
and  standard  deviation.  These  features  are  similar  for  the  model  errors  on  the  180th  day. 
These  relatively  strong  patterns  are  also  present  in  Run  11,  even  though  the  difference 
between  Run  1 1  and  9  is  not  quite  obvious.  This  indicates  the  model  errors  caused  by 
Run  1 1  have  a  larger  effect  on  the  model  performance  than  that  caused  by  Run  9  and  10. 
As  observed,  the  model  errors  are  present  mostly  near  the  Japan  Basin.  This  could  be 
associated  with  the  nearby  Polar  Front  (PF)  where  two  current  systems  (NKCC  and 
TWC)  converge.  In  general,  from  these  horizontal  distribution  plots,  obvious  differences 
exist  between  each  experimental  run  (the  maximum  and  minimum  model  error  is  caused 
by  Run  1 1  and  Run  10,  respectively)  with  the  model  errors  decreasing  obviously  from  the 
5th  to  the  180th  day. 
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b.  The  Histogram 

The  model  error,  calculated  by  equation  (11),  could  be  presented  in  the 
form  of  a  histogram  in  evaluating  the  model  error  distribution.  Three  histogram  plots 
contained  the  model  errors  on  each  specific  field  (u  and  v)  between  the  control  run  and 
three  experimental  runs  made  on  each  5th  day  and  180th  day.  The  values  of  the  number  of 
the  total  sample,  samples  greater  than  zero,  samples  less  than  zero  and  standard  deviation 
were  also  determined  for  each  plot.  These  plots  are  shown  in  Figure  32  to  33. 

In  Figure  32,  the  histogram  of  model  errors  of  u  components  of  the 
horizontal  velocity  on  both  the  5th  day  and  the  180th  day  also  reveals  a  Gaussian-type 
distribution  on  each  plot.  Similar  to  previous  plots,  the  values  of  standard  deviation 
(STD)  are  higher  in  Run  11  than  those  in  Run  9  and  10,  and  they  decrease  with 
integration  time  on  each  experimental  run.  For  the  5  th  day,  the  maximum  and  minimum 
values  of  standard  deviation  occur  in  Run  1 1  and  Run  10  with  a  value  of  0.026  m/s  and 
0.022  m/s,  respectively.  For  the  180th  day,  the  maximum  and  minimum  values  of  standard 
deviation  also  occur  in  Run  11  and  Run  10  with  a  value  of  0.018  m/s  and  0.01  m/s, 
respectively.  These  features  are  similar  for  the  histogram  of  model  errors  of  v 
components  of  the  horizontal  velocity  on  both  the  5th  day  and  the  180th  day  in  Figure  33. 
Once  again,  for  both  the  5th  and  the  180th  day,  the  maximum  and  minimum  values  of 
standard  deviation  occur  in  Run  11  and  Run  10  with  the  maximum  and  minimum  value 
of  0.029  m/s  and  0.024  m/s  on  the  5th  day  and  0.02  m/s  and  0.012  m/s  on  the  180th  day, 
respectively.  General  speaking,  the  model  error  caused  by  Run  1 1  has  the  largest  effect 
on  model  performance  among  these  three  experimental  runs. 
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Figure  30.  Vertically  averaged  model  error  of  (a)  Run  0  -  Run  9,  (b)  Run  0  -  Run  10 
and  (c)  Run  0  -  Run  11  for  u  component  of  velocity  field  on  both  the  5th 
and  180th  day  after  the  model  integration. 


60 


Latitude  Latitude  Latitude 


180th  Day 


5th  Day 


-0.05  -0.04  -0.03  -0.02  -0.01  0  0.01  0.02  0.03  0.04  0.05 

Figure  3 1 .  Vertically  averaged  model  error  of  (a)  Run  0  -  Run  9,  (b)  Run  0  -  Run  10 
and  (c)  Run  0  -  Run  11  for  v  component  of  velocity  field  on  both  the  5th 
and  180th  day  after  the  model  integration. 
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Figure  32.  Histogram  of  (a)  Run  0  -  Run  9,  (b)  Run  0  -  Run  10  and  (c)  Run  0  -  Run 
11  for  u  component  of  velocity  field  on  both  the  5th  and  180th  day  after  the 
model  integration. 
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Figure  33.  Histogram  of  (a)  Run  0  -  Run  9,  (b)  Run  0  -  Run  10  and  (c)  Run  0  -  Run 
11  for  v  component  of  velocity  field  on  both  the  5th  and  180th  day  after  the 
model  integration. 
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2. 


The  Root  Mean  Square  Error  (RMSE) 


a.  The  Vertical  Variation 

The  Root  Mean  Square  Error,  calculated  by  equation  (12),  is  used  to 
evaluate  the  model  perfonnance.  For  the  vertical  profile  of  RMSE,  two  plots  contained 
the  errors  caused  by  these  three  experimental  runs  for  the  horizontal  velocity  (V)  fields 
made  on  both  the  5th  day  (day- 185)  and  the  180th  day  (day-360)  after  the  day  of  model 
integration  (day- 180).  These  plots  are  shown  in  Figure  34. 

As  seen  in  this  figure,  the  RMSE  on  the  5th  day  reveals  that  a  higher  value 
occurs  above  sigma  level  8,  which  indicates  that  some  stronger  effects  on  model 
prediction  occur  at  those  levels.  The  maximum  RMSE  occurs  in  Run  1 1  at  level  1  with  a 
value  of  0.08  m/s.  These  features  are  similar  for  the  RMSE  on  the  180th  day.  The 
maximum  RMSE  occurs  in  Run  1 1  at  level  1  with  a  value  of  0.075  m/s.  Although  no 
obvious  differences  take  place  between  Run  9  and  11,  the  RMSE  caused  by  Run  1 1  has 
the  largest  value  among  these  three  experimental  runs.  In  addition,  these  values  of  RMSE 
decrease  with  depth  (rapidly  above  sigma  level  8). 

b.  The  Temporal  Evolution 

The  temporal  evolution  of  vertically  averaged  RMSE  for  horizontal 
velocity  (V)  fields  from  the  5th  day  to  the  180th  day  is  shown  in  Figure  35.  As  seen  in  this 
figure,  a  significant  difference  exists  among  these  three  curves  of  RMSE  for  the 
horizontal  velocity  (V)  fields  (more  RMSE  is  caused  by  Run  11).  As  observed,  the 
RMSE  is  largely  detennined  due  to  the  initial  uncertainty  in  the  early  period  (before  the 
45th  day),  the  lateral  boundary  transport  uncertainty  in  the  middle  period  (between  the 
45th  day  and  the  155th  day)  and  the  uncertain  winds  in  the  late  period  (near  the  180th  day). 
The  RMSE  of  horizontal  velocity  oscillates  and  generally  decreases  with  time  with  a 
peak  value  on  the  5th  day  of  0.035,  0.034  and  0.031  m/s  in  Run  11,  9  and  10,  respectively. 
As  in  previous  results,  the  higher  RMSE  is  caused  by  Run  1 1  than  by  Run  9  and  10. 
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Figure  34.  RMSE  of  V  in  (a)  Run  0  -  Run  9,  (b)  Run  0  -  Run  10  and  (c)  Run  0  -  Run 
1 1  on  the  5th  day  and  the  180th  day  after  the  model  integration. 


Figure  35.  Temporal  evolution  of  vertically  averaged  RMSE  of  V  in  (a)  Run  0  -  Run 
9,  (b)  Run  0  -  Run  10  and  (c)  Run  0  -  Run  1 1 . 
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3. 


The  Relative  Root  Mean  Square  Error  (RRMSE) 


a.  The  Vertical  Variation 

The  Relative  Root  Mean  Square  Error,  calculated  by  equation  (13),  is  also 
commonly  used  for  evaluating  the  model  performance.  For  the  vertical  profile  of 
RRMSE,  two  plots  contained  the  errors  caused  by  these  three  experimental  runs  for  the 
horizontal  velocity  (V)  fields  made  on  both  the  5th  day  (day- 185)  and  the  180th  day  (day- 
360)  after  the  day  of  model  integration  (day-180).  These  plots  are  shown  in  Figure  36. 

As  seen  from  this  figure,  the  RRMSE  on  the  5th  day  reveals  that  the 
relatively  higher  RRMSE  occurs  in  Run  1 1 .  The  maximum  RRMSE  occurs  in  Run  1 1  at 
level  1  with  a  value  of  0.73.  For  the  RRMSE  on  the  180th  day,  it  also  reveals  that  the 
relatively  higher  RRMSE  of  all  three  fields  occurs  in  Run  11.  Except  for  Run  10,  a  higher 
value  of  RRMSE  is  present  above  sigma  level  8,  which  indicates  some  stronger  effects  on 
model  prediction  occurring  at  those  levels.  The  maximum  RRMSE  occurs  in  Run  1 1  with 
a  value  of  0.78  at  level  1  for  horizontal  velocity.  As  in  previous  plots,  although  no 
obvious  differences  exist  between  Run  9  and  1 1 ,  the  RRMSE  caused  by  Run  1 1  has  a 
relatively  higher  value  among  these  three  experimental  runs.  Generally,  the  RRMSE 
decreases  with  depth  (rapidly  above  sigma  level  8)  except  for  Run  10.  As  observed,  the 
RRMSE  of  horizontal  velocity  in  Run  1 1  is  quite  significant,  especially  above  sigma 
level  8  (around  73  %  and  78  %  on  the  5th  day  and  the  180th  day  at  the  surface  level, 
respectively). 


b.  The  Temporal  Evolution 

The  temporal  evolution  of  vertically  averaging  RRMSE  for  horizontal 
velocity  (V)  fields  from  the  5th  day  to  the  180th  day  is  shown  in  Figure  37.  As  seen  in  this 
figure,  significant  differences  exist  among  these  three  curves  of  RRMSE  for  the 
horizontal  velocity  (V)  fields  (more  RRMSE  was  caused  by  Run  11).  As  observed,  the 
RRMSE  is  largely  detennined  due  to  the  initial  uncertainty  in  the  early  period  (before  the 
40th  day),  the  lateral  boundary  transport  uncertainty  in  the  middle  period  (between  the 
40th  day  and  the  160th  day)  and  the  uncertain  winds  in  the  late  period  (near  the  180th  day). 
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The  RRMSE  of  horizontal  velocity  oscillates  and  generally  decreases  with  time  with  a 
peak  value  on  the  5th  day  of  0.52,  0.5  and  0.48  in  Run  11,9  and  10,  respectively.  As  in  all 
previous  results,  the  higher  RRMSE  is  caused  by  Run  1 1  than  by  the  Run  9  and  10. 
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Figure  36.  RRMSE  of  V  in  (a)  Run  0  -  Run  9,  (b)  Run  0  -  Run  10  and  (c)  Run  0  - 
Run  lion  the  5  th  day  and  the  180th  day  after  the  model  integration. 


Figure  37.  Temporal  evolution  of  vertically  averaged  RRMSE  of  V  in  (a)  Run  0  - 
Run  9,  (b)  Run  0  -  Run  10  and  (c)  Run  0  -  Run  1 1 . 
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VII.  CONCLUSIONS 


In  this  thesis  several  experiments  were  conducted  in  order  to  investigate  the 
model  errors  due  to  input  uncertainty  of  the  velocity  initial  conditions,  wind  forcing  and 
lateral  boundary  transport.  The  RRMSE  of  the  horizontal  velocity  field  in  these 
experiments  is  summarized  in  Table  7. 


Experiment 

Vertically  averaged 
RRMSE 

Max.  RRMSE 

Min. 

Max. 

5th  Day 

180th  Day 

For  uncertain  velocity 
initial  conditions 

20% 

50% 

70%  near 
the  surface 

25%  near 
the  surface 

For  0.5  m/s  noise 
intensity 

8% 

19% 

35%  near 
the  surface 

50%  near 
the  surface 

For  1.0  m/s  noise 
intensity 

11% 

28% 

60%  near 
the  surface 

80%  near 
the  surface 

For  noise  intensity  as 

5%  of  transport 

9% 

20% 

14%  near 
the  bottom 

18%  near 
the  bottom 

For  noise  intensity  as 

10%  of  transport 

17% 

34% 

24%  near 
the  bottom 

28%  near 
the  bottom 

For  uncertain  initial 
condition  and  wind 
forcing 

20% 

52% 

70%  near 
the  surface 

77%  near 
the  surface 

For  uncertain  initial 
condition  and  lateral 
boundary  transport 

27% 

50% 

65%  near 
the  bottom 

35%  near 
the  bottom 

For  uncertain  initial 
condition,  wind  forcing 
and  lateral  boundary 
transport 

30% 

55% 

73%  near 
the  surface 

78%  near 
the  surface 

Table  7  Summary  of  the  RRMSE  in  each  experiment. 
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For  uncertain  velocity  initial  conditions  with  and  without  diagnostic  initialization, 
the  model  errors  are  quite  comparable  and  significant,  but  they  decrease  with  time.  For 
diagnostic  initialization,  the  magnitude  of  model  errors  is  less  dependent  on  the 
initialization  period  no  matter  if  it  is  30  days  or  90  days.  The  vertically  averaged  RRMSE 
of  the  horizontal  velocity  decreases  rapidly  from  50%  (the  maximum  value)  on  the  5th 
day  to  20%  (the  minimum  value)  on  the  50th  day  and  then  slightly  oscillates  with  time 
near  20%  to  the  180th  day.  The  RRMSE  of  horizontal  velocity  is  70%  near  the  surface 
and  50%  near  the  bottom  on  the  5th  day  and  25%  near  the  surface  and  20%  near  the 
bottom  on  the  180th  day. 

For  uncertain  wind  forcing  with  the  Gaussian  random  noise,  the  model  error 
increases  with  time  and  noise  intensity.  The  vertically  averaged  RRMSE  of  the  horizontal 
velocity  fluctuates  with  time.  For  the  noise  intensity  of  0.5  m/s,  it  increases  slowly  from 
1 1%  on  the  5th  day  to  12%  on  the  45th  day  and  then  decreases  to  8%  (the  minimum  value) 
on  the  75th  day  and  then  increases  again  to  19%  (the  maximum  value)  on  the  180th  day. 
For  the  noise  intensity  of  1.0  m/s,  it  increases  slowly  from  18%  on  the  5th  day  to  20%  on 
the  45th  day  and  then  decreases  rapidly  to  11%  (the  minimum  value)  on  the  75th  day  and 
then  increases  again  to  28%  (the  maximum  value)  on  the  180th  day.  The  maximum 
RRMSE  of  horizontal  velocity,  occurring  near  the  surface,  increases  from  35%  on  the  5th 
day  to  50%  on  the  180th  day  for  noise  intensity  of  0.5  m/s,  and  increases  from  60%  on  the 
5th  day  to  80%  on  the  180th  day  for  noise  intensity  of  1.0  m/s.  The  model  errors  generally 
decrease  with  depth  and  become  small  below  sigma  level-8. 

For  uncertain  lateral  boundary  transport  with  the  Gaussian  random  noise,  the 
model  error  increases  with  time  and  noise  intensity.  The  vertically  averaged  RRMSE  of 
the  horizontal  velocity  fluctuates  with  time  between  9%  and  20%  (17%  and  34%)  for  the 
noise  intensity  to  be  5%  (10%)  of  the  transport  of  the  control  run.  It  generally  increases 
with  time  from  the  5th  day  to  the  150th  day  with  a  peak  value  of  20%  for  5%  noise  and 
34%  for  10%  noise  and  then  decreases  with  time  to  15%  for  5%  noise  and  to  22%  for 
10%  noise  on  the  180th  day.  The  maximum  RRMSE  of  horizontal  velocity,  occurring 
near  the  bottom,  increases  from  14%  on  the  5th  day  to  18%  on  the  180th  day  for  5%  noise, 
and  increases  from  24%  on  the  5th  day  to  28%  on  the  180th  day  for  10%  noise.  The  model 
errors  generally  increase  with  depth. 
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For  combined  uncertain  initial  condition  (30  day  period  diagnostic  initialization) 
and  wind  forcing  (1.0  m/s  noise  intensity),  the  vertically  averaged  RRMSE  of  the 
horizontal  velocity  fluctuates  between  52%  and  20%  while  decreasing  from  52%  on  the 
5th  day  to  20%  on  the  50th  day  and  then  increasing  slowly  to  32%  on  the  180th  day.  The 
model  uncertainty  is  largely  determined  due  to  the  initial  uncertainty  in  the  early  period 
(before  the  50th  day)  and  due  to  the  uncertain  winds  in  the  late  period  (near  the  180th  day). 
The  maximum  RRMSE  of  horizontal  velocity,  occurring  near  the  surface,  increases  from 
70%  on  the  5th  day  to  77%  on  the  180th  day. 

For  combined  uncertain  initial  condition  (30  day  period  diagnostic  initialization) 
and  lateral  boundary  transport  (10%  noise  intensity),  the  vertically  averaged  RRMSE  of 
the  horizontal  velocity  fluctuates  between  50%  and  27%  while  decreasing  from  50%  on 
the  5th  day  to  30%  on  the  30th  day  and  then  oscillates  between  27%  and  38%  with  time  to 
29%  on  the  180th  day.  The  model  uncertainty  is  largely  determined  due  to  the  initial 
uncertainty  in  the  early  period  (before  the  50th  day)  and  the  lateral  boundary  transport 
uncertainty  in  the  late  period  (near  the  150th  day).  The  maximum  RRMSE  of  horizontal 
velocity,  occurring  near  the  bottom,  decreases  from  65%  on  the  5th  day  to  35%  on  the 
180th  day. 

For  combined  uncertain  initial  condition  (30  day  period  diagnostic  initialization), 
wind  forcing  (1.0  m/s  noise  intensity)  and  lateral  boundary  transport  (10%  noise 
intensity),  the  vertically  averaged  RRMSE  of  the  horizontal  velocity  fluctuates  between 
55%  and  30%  while  decreasing  from  55%  on  the  5th  day  to  36%  on  the  30th  day  and  then 
oscillates  between  30%  and  45%  with  time  to  38%  on  the  180th  day.  The  model 
uncertainty  is  largely  determined  due  to  the  initial  uncertainty  in  the  early  period  (before 
the  50th  day),  the  lateral  boundary  transport  uncertainty  in  the  middle  period  (between  the 
50th  day  and  the  150th  day)  and  the  uncertain  winds  in  the  late  period  (near  the  180th  day). 
The  maximum  RRMSE  of  horizontal  velocity,  occurring  near  the  surface,  increases  from 
73%  on  the  5th  day  to  78%  on  the  180th  day. 

As  the  table  demonstrates,  the  RRMSE  with  such  higher  values  in  each 
experiment  reveal  that  the  model  uncertainty  due  to  these  input  uncertainties  (the  velocity 
initial  conditions,  wind  forcing  and  lateral  boundary  transport)  is  very  significant. 
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Therefore,  it  is  necessary  to  reduce  these  uncertainties  to  achieve  better  prediction  of  the 
ocean  behavior.  In  addition,  the  results  in  this  thesis  can  only  be  applied  to  this  specific 
numerical  ocean  model  (POM),  study  area  (JES)  and  model  setup  (i.e.,  resolution). 
Applying  these  results  to  other  ocean  models,  areas  or  model  setups  needs  further 
investigation. 
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